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Abstract 

Dark-matter (DM) halos are the scaffolding around which galaxies and clus- 
ters are built. They form when the gravitational instability of primordial den- 
sity fluctuations causes regions which are denser than average to slow their 
cosmic expansion, recollapse, and virialize. Understanding the equilibrium 
structure of these halos is thus a prerequisite for understanding galaxy and 
cluster formation. Numerical N-body simulations of structure formation from 
Gaussian-random-noise initial conditions in the Cold Dark Matter (CDM) uni- 
verse find a universal internal structure for halos. Objects as different in size 
and mass as dwarf spheroidal galaxies and galaxy clusters are predicted to have 
halos with the same basic structure when properly rescaled, independent of halo 
mass, of the shape of the power spectrum of primordial density fluctuations, 
and of the cosmological background parameters. This remarkable universality 
is a fundamental prediction of the CDM model, but our knowledge is limited to 
the "empirical" N-body simulation results, with little analytical understanding. 
We summarize here our attempts to fill this gap, in an effort to derive and give 
physical insight to the numerical results and extend them beyond the range of 
numerical simulation: (1) Simulated halos which form from highly simplified 
initial conditions involving gravitational instability in a cosmological pancake 
show that many of the universal properties of CDM halos are generic to cosmo- 
logical gravitational collapse and do not require Gaussian-random-noise density 
fluctuations or hierarchical clustering. (2) A fluid approximation derived from 
the Boltzmann equation yields an analytical theory of halo dynamics which can 
explain many of the N-body results if the complex mass assembly history of indi- 
vidual halos is approximated by continuous spherical infall. The universal mass 
growth history reported for CDM N-body halos corresponds to a time-varying 
infall rate which self-consistently determines the shape of the equilibrium halo 
profile and its evolution without regard for the complicated details of the merger 
process. (3) The first fully-cosmological, similarity solutions for halo formation 
in the presence of collisionality provide an analytical theory of the effect of the 
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self-interacting dark matter (SIDM) hypothesis on CDM halo density profiles. 
Collisions transport heat inward which flattens the central cusp, but continu- 
ous infall pumps energy into the halo to stabilize the core against gravother- 
mal catastrophe. (4) The postcollapse equilibrium halo structure derived by 
matching spherical top-hat collapse to a unique, minimum-energy solution of 
the isothermal Lanc-Emdcn equation the Truncated Isothermal Sphere (TIS) 
model - reproduces well the average properties of simulated CDM halos. The 
TIS yields a universal, self-similar equilibrium structure whose parameters are 
determined by the halo's total mass and collapse redshift. These latter two 
parameters are statistically correlated, however, since halos of the same mass 
form on average at the same epoch, with small-mass objects forming first and 
then merging hierarchically to form larger-mass objects. The structural prop- 
erties (e.g. sizes, mass profiles, velocity dispersions and virial temperatures) 
of dark-matter dominated halos of different masses, therefore, should reflect 
this statistical correlation, an imprint of the statistical properties of the pri- 
mordial density fluctuations which formed them. When combined with the 
Press-Schcchter (PS) formalism, the TIS model can predict these statistical 
correlations of halo properties analytically, which we compare to observational 
data on these correlations, providing a fundamental test of the CDM model, 
which probes the shape of the power spectrum of primordial density fluctua- 
tions and the cosmological background parameters. (5) For observational tests 
which are sensitive to the shape of the inner mass profile at the very center, 
the TIS model provides a convenient analytical tool for studying the effect of a 
small nonsingular core on CDM halo predictions. As an example, strong gravi- 
tational lensing properties of CDM halos with nonsingular cores are contrasted 
with those with singular, cuspy cores. 



1 Introduction 

Numerical N-body simulations of structure formation in the Cold Dark Matter 
(CDM) universe show a cosmic web of lumps, filaments, and sheets. This web 
evolves in a self-similar way, with smaller mass structures forming first and merging 
to form larger-mass structures later, in a continuous sequence of mass assembly. 
At each epoch, the web produces gravitationally-bound, quasi-spherical "halos" in 
virial and hydrostatic equilibrium. These virialized halos are the sites of galaxy and 
cluster formation. Their universal equilibrium structure is a fundamental prediction 
of the CDM model, but our knowledge is limited to the numerical N-body results, 
with very little analytical understanding. We shall describe our attempts to fill this 
gap in what follows. 

Most of the progress to date on the formation and evolution of virialized, dark- 
matter dominated halos in a CDM universe has been via numerical N-body simula- 
tions of collisionless dark matter involving Gaussian random noise initial conditions. 
According to those N-body simulations, the spherically-averaged mass distribution 
inside halos is universal, with a density profile which declines with radius, approach- 
ing p oc r -3 at large radii, flattening near the center to p(r) oc r~ a with a < 2. Two 
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"universal" profiles bracket the results (NFW, QUE]; Moore, [100]). The NFW 
(Moore) profile has an inner density profile p oc r (r ). The actual value of 
a is still uncertain, even though the N-body results have advanced to the point of 
including millions of simulation particles within the virial radius of a given halo (e.g. 
[54* 1 1551 I5H 1 158*1 ITil I75 | "8l"l"ll2| ). Since the true inner profile may not be a power-law, 
comparisons of the inner slope for different simulation results are typically referred 
to a particular radius (e.g. r = 0.01r v ; r , where r v - u is some measure of the outer 
radius of the virialized region). The N-body results generally support the conclusion 
that the same halo density profile applies to objects as different as dwarf galaxies 
and galaxy clusters, independent of halo mass, of the shape of the density fluctuation 
power spectrum, and of the background cosmology. This universality is apparent 
when comparisons of different halo profiles are made with density expressed in units 
of the density p_2 = p{f—2) at the radius r_2 at which each profile has a logarith- 
mic slope of —2, while radii for each halo are expressed in units of r_2 109 . An 
exception to this universality is claimed by |119j . who reports that the value of a 
depends upon halo mass, shallower for dwarf galaxies than for clusters, reflecting the 
different slopes of the power spectrum at the different scales represented by these 
objects, but [HH] has challenged this claim. 

Much attention has been focused on the N-body results for this inner slope, 
since the observed rotation curves of dark-matter dominated dwarf and low surface 
brightness (LSB) disk galaxies tend to favor mass profiles with a flat-density core 
unlike the singular profiles of the CDM N-body simulations (e.g. [Tfl t88l I5TI t98l fT3Tj : 
but, for a different view, see also [138, 139]). On the cluster scale, too, there is some 
evidence from observations of strong gravitational lensing of background galaxies by 
foreground clusters which favors a flatter inner halo density profile than is found by 
the CDM N-body simulations (e.g. "57] Q231 ESI ; but see also EH)- However, the 
halo mass fraction contained within the disputed inner cusp is actually quite small, 
so there may yet be dynamical processes not fully accounted for in the pure N-body 
simulations which can affect this small central mass without disturbing the overall 
universality of the rest of the profile. 

Along with their universal mass profiles, CDM N-body halos also exhibit several 
universal properties in their phase-space distributions. Over most of the halo volume 
inside the virial radius, the DM particles are approximately isothermal - i.e. their 
velocity dispersion is fairly constant with radius - with only a relatively small dip 
in "temperature" toward the center (e.g. |52t 1144] ). Halo particle velocities are also 
approximately isotropic, with only a mild radial bias in the outer halo, which gives 
way to increasing isotropy toward the center (e.g. [2I1I"""DII"3[5I]). The spherically- 
averaged mass motion at each radius is quite small and "subsonic;" the halo is not 
only in a state of global virial equilibrium but is close to hydrostatic equilibrium at 
each radius, too - i.e. it satisfies the spherical Jeans equation (e.g. |144j ). 

Individual halos in CDM N-body simulations evolve over time, on average, 
through a continuous sequence of universal-shaped mass profiles of increasing to- 
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tal mass [14=11 11491 I15U) . This Lagrangian mass evolution can be characterized by 
a universal mass accretion history: M(a)/M(a,f), where a is the cosmic scale factor 
and cif is some particular value of a, such as that at which din Mj dm a = 2 150 . 
As the mass of each halo grows with time due to the average effect of mergers and 
smooth infall, so does the concentration parameter c of its density profile, where 
c = r v ; r /r_2> roughly as c(a)/c(a,f) oc a/at for a/af > 1 j!5Uj . after hovering at low 
values c < 3 — 4 during the initial phase of most rapid mass assembly prior to at 

cm- 

This description of the CDM halos of N-body simulations is a spherically-averaged 
one, which neglects many details. There is some scatter in the N-body results about 
this average description, of course. Individual halos are not truly spherically sym- 
metric, either, but only approximately so. The neglect of net angular momentum 
is probably not a bad first approximation, since the specific angular momentum is 
typically found to be far below that required for rotational support (e.g. [121 .22 J 1 . 
However, the spherically- averaged description also averages out the small-scale den- 
sity inhomogeneities inside each halo. This small-scale inhomogeneity may play an 
important role in the underlying dynamics which leads to halo formation and evolu- 
tion in these N-body simulations. Regardless of its dynamical significance, this halo 
substructure has also been the subject of special attention for its own sake, once 
it was noticed that the number of subhalos which typically survive their merger 
into a larger halo in the N-body results is much larger than the number of galaxies 
observed within the Local Group (e.g. [80" 1 1101] ) . 

While N-body simulations have made the universal equilibrium structure of ha- 
los described above a fundamental prediction of the CDM model, much less progress 
has been made on the analytical side, to derive and understand the numerical results 
and extend them beyond the range of numerical simulation. Our purpose in what 
follows is to summarize our own attacks on this problem by a hierarchy of approxi- 
mations, each one simpler than the last, from 3D gas and N-body simulations of halo 
formation involving simplified initial conditions (©, to ID, spherically-symmetric, 
dynamical models involving a fluid approximation derived from the Boltzmann equa- 
tion (^HJ), to a model for the hydrostatic equilibrium of halos which follows from the 
virialization of top-hat density perturbations (HU). As we shall see, the last of these, 
known as the Truncated Isothermal Sphere (TIS) model, not only provides an ex- 
cellent match to many of the average properties of the halos found by the N-body 
simulations of CDM from realistic initial conditions, but one which is conveniently 
coupled to the Press-Schechter formalism for the mass function of halos in CDM 
to predict the observed statistical properties of galaxies and clusters analytically. 
We summarize several such applications of this model in ^4.31 In motivated by 
the prospect that gravitational lensing will provide a direct observational test of 

lr The specific angular momentum profile, j{M), of individual CDM halos averaged over spherical 
shells encompassing mass M, also has a universal shape, which has been fitted by j(M) oc M" with 

s = i.3 ±0.3 Eg. 
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the CDM halo density profile, we use the analytical theory of the TIS halo, with 
its small, flat-density core at the center of a halo which otherwise closely resembles 
the average CDM halos of N-body simulations, to contrast the lensing properties of 
nonsingular and singular CDM halos. 

2 N-body Halos from Simplified Initial Conditions: 3D 
Halo Formation by Gravitational Instability of Cos- 
mological Pancakes 

A natural question which emerges from the results of N-body simulations of struc- 
ture formation in CDM is whether the universal properties of dark matter halos, 
like their density profiles and mass accretion histories, are a consequence of hier- 
archical clustering from Gaussian-random-noise density fluctuations or are in fact 
more general. To address this question, we have analyzed the results of simula- 
tions with initial conditions that are much simpler than a CDM power spectrum 
with Gaussian-random-noise, while retaining the realistic features of continuous and 
anisotropic infall in three dimensions. Each simulation forms a single halo from the 
gravitational instability of a cosmological pancake. Because we focus on one halo at 
a time, we are able to follow its formation and evolution, and find that several of the 
trends reported for halo evolution in the CDM simulations are also present in the 
simplified pancake model. Thus, the pancake instability model serves as a conve- 
nient test-bed for halo formation which shares not only similar boundary conditions 
with more realistic models of halo formation, but also produces halos that have 
similar structural and evolutionary properties. Some of this work was summarized 
previously in jSJ H3 Ej ■ 

2.1 Halo Formation via Pancake Instability 

2.1.1 Unperturbed Pancake 

Consider the growing mode of a single sinusoidal plane-wave density fluctuation of 
comoving wavelength A p and dimensionless wavevector k p = x (length unit = A p ) 
in an Einstein-de Sitter universe dominated by cold, collisionless dark matter |13flj . 
Let the initial amplitude 6\ at scale factor a\ be chosen so that the first density 
caustic forms in the collisionless component at scale factor a = a c = a{/5\. 

2.1.2 Perturbations 

Pancakes modeled in this way have been shown to be gravitationally unstable, lead- 
ing to filamentation and fragmentation during the collapse |148| . As an example, we 
shall perturb the ID fluctuation described above by adding to the initial primary 
pancake mode two transverse, plane- wave density fluctuations with equal wavelength 
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Figure 1: Dark matter particles at o/a c = 3. (see jpg file) 



A s = A p , wavevectors k s pointing along the orthogonal unit vectors y and z, and 
smaller initial amplitudes, e y 8\ and e z 8\, respectively, where e 9 «l and f 2 « 1. A 
pancake perturbed by two such density modes will be referred to as Sx >eyj€z . All re- 
sults presented here refer to the case 51,0.2,0.2 ■ The initial position, velocity, density, 
and gravitational potential are given by 

n = qi + ^jjr sin ^k p q t i 1 ) 

ei ( d5\ 



and 



= P , ( 3 ) 

1 + 5i(cos 2irk p q x + e y cos 2nk p q y + e z cos 2-7rk p q z ) ' 

4> = (4>) (cos 2irk p x + e y cos 2nk p y + e z cos 2i:k p z) , (4) 



where I = x, y, or z, e x = 1, q x , q y , and q z are the unperturbed particle positions, 
and ~p is the cosmic mean matter density. 

Such a perturbation leads to the formation of a quasi-spherical mass concentra- 
tion in the pancake plane at the intersection of two filaments (Fig. 0). 
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2.2 N-body and Hydrodynamical Simulations 

The code we use to simulate the formation of the halo couples the Adaptive SPH 
(ASPH) algorithm, first described in Q29 , to a P 3 M gravity solver. The 
ASPH method improves on standard SPH by introducing nonspherical, ellipsoidal 
smoothing kernels to better track the anisotropic flow that generally arises during 
cosmological structure formation. Two simulations were carried out, one with gas 
and one without. In both cases, there were 64 3 particles of dark matter, while 
there were also 64 3 gas particles when gas was included. The P 3 M grid was 128 
cells in a periodic cube size A p on a side, with a comoving softening length of 
r soft = 0.3Ax = 0.3A p /128, where Ax is the cell size. The initial conditions were 
those described in § 12.1.21 

The adiabatic pancake problem (i.e. without radiative cooling) is self-similar 
and scale-free, once distance is expressed in units of the pancake wavelength A p and 
time is expressed in terms of the cosmic scale factor a in units of the scale factor 
a c at which caustics first form in the dark matter and shocks in the gas ^30 In 
the currently-favored, flat, cosmological-constant-dominated universe, however, this 
self-similarity is broken because Om/^a decreases with time, where Om and Oa are 
the matter and vacuum energy density parameters, respectively. For objects which 
collapse at high redshift in such a universe (e.g. dwarf galaxies), the Einstein-de 
Sitter results are still applicable as long as we take (Ob/Odm)e<is = (^b/^dm)A) 
where Ob and Odm ar e the baryon and dark matter density parameters. If Ob = 
0.045, Odm = 0.255, and Oa = 0.7 at present, then the EdS results are applicable 
if we take Ob = 0.15 and Odm = 0.85, instead. 



2.3 Profiles 
2.3.1 Density 

The density profiles at different epochs for the simulations without gas are shown in 
Fig. along with the best-fitting NFW profile for each epoch, which has the form 

" * (5) 



with 5 C given by 



p (r/r s )(l + r/r s ) 2 ' 
200 c 3 



3 ln(l + c) -c/(l + c)' (6) 

where r s = r_2, c = r2oo/r s , and r2oo is the radius within which the mean density 
(p) = 200/5. The NFW profile has only one free parameter, the concentration pa- 
rameter c. Our pancake halo density profiles have concentrations which range from 
3 to 15, increasing with time, and are usually within 20% of the best-fit NFW profile 
at all radii. This NFW profile is a fit to N-body results for CDM halos, but it is also 
consistent with the halos which form in simulations using Gaussian-random-noise 




Figure 2: Density profile of the dark matter halo as simulated without gas at four 
different scale factors, a/a c =3, 5, 7, and 10, with spherically-averaged simulation 
results in radial bins (filled circles) and the best-fitting NFW profiles (solid curves) 
for several epochs, as labeled. Shown above each panel are fractional deviations 
(/Onfw — 9)1 Pnfw from the best-fitting NFW profiles for each epoch. Vertical lines 
indicate the location of r so f t , the numerical softening-length, and r2oo- 



initial fluctuations with small-scale fluctuations suppressed (HI E21 EUl E2] , as is 
the case for the halos that form from pancake instability presented here. 
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Figure 3: Dimensionless specific thermal energy profiles at a/a c = 7 in the 
(gas+DM)-simulation, for the dark matter (left) and gas (right), with dotted lines 
indicating the rms scatter within each logarithmic bin. 

2.3.2 Kinetic and Thermal Energy 

The profiles of the dimensionless specific thermal kinetic energy 6dm = (l v — 
(v)| 2 )/2 = 3<7£) M /2 of the dark matter and the dimensionless specific thermal energy 
e gas = 3k^T/2m of the gas are shown in Fig. 01 f° r the simulation with gas and dark 
matter at a/a c = 7. Although the dark matter velocity dispersion rises towards 
the center, the rise is shallow and the kinetic energy distribution is approximately 
isothermal inside the radius r2oo- As seen from Fig. El the gas is even more isother- 
mal than the dark matter, with e gas (r so f t )/e gas (r2oo) — 2.5, while the density varies 
by more than three orders of magnitude over the same region. In Fig. 0]we show 
the ratio of the specific thermal energy of the dark matter to that of the gas. This 
ranges from eoMAgas ~ 1-6 in the center to eDiiAgas ~ 1 at r2oo- 

2.3.3 Velocity Anisotropy 

In Fig. 01 we plot the profile of the anisotropy parameter 0, defined two different 
ways, according to the frame of reference in which the velocity dispersion is calcu- 
lated. In the Eulerian case, where the bulk motion of the shell contributes to the 
anisotropy, it is defined as = 0e = 1 — ( v t) / (^( v r)) ■ I* 1 the Lagrangian case, 
however, the bulk motion of the shell is subtracted out, = 0i, = 1 — of/(2cr 2 ), 
where of = ({vi — (vi)) 2 ). Because of the symmetry of the pancake initial con- 
ditions, (v 2 ) = of at all times. Differences between the two definitions therefore 
arise because of differences in the radial component. The values = 1,0, and — oo 
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Figure 4: (left) Ratio of dimensionless specific thermal energy profiles at a/a c = 7 
in the (gas+DM)-simulation. (right) Anisotropy profile, defined two different ways, 
at a/o c = 7. The curves are labeled according to the definitions given in the text. 



correspond to motion which is purely radial, fully isotropic, and tangential, respec- 
tively. While the two definitions of anisotropy give nearly indistinguishable profiles 
at r < r200; the two profiles depart significantly at r > V2oq- This is expected, 
since the bulk radial motion is nearly zero inside the halo where equilibrium is a 
reasonable expectation. The two definitions of (3 are identical when the system is 
in equilibrium. Outside the halo, however, equilibrium is violated because the bulk 
radial motion is not zero, and the difference between the profiles arises due to the 
detailed nature of the region outside the halo. 

Simulations of CDM typically find values of (3 near at the center, slowly rising 
to a value of (3 « 0.5 at r2oo EZ1 EU ESI EH- As seen from Fig. |1J the halo formed 
by pancake instability is more anisotropic, with values of (3 rising from ~ 0.2 near 
the center to ~ 0.8 at ^200- This reflects the strongly filamentary substructure of 
the pancake within which the halo forms, and perhaps the absence of strong tidal 
fields or mergers as well, which might otherwise help convert radial motions into 
tangential ones. 

2.4 Virial Equilibrium 

A state of equilibrium is commonly assumed in analytical modeling of dark matter 
halos [22] ■ Using A-body simulations, |144j found that halos which formed from 
CDM initial conditions roughly obey the Jeans equation for dynamical equilibrium 
in spherical symmetry, within a radius of order r2oo> suggesting that CDM halos are 
in approximate virial equilibrium. In what follows, we will interpret the numerical 
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Figure 5: (left) Virial ratio versus dimensionless radius ( for the TIS solution. The 
vertical dotted line on the right corresponds to the truncation radius, Q, at which 
the total energy is a minimum at fixed mass and boundary pressure, while the one 
on the left corresponds to the radius within which the mean density is 200 times 
the background density. The horizontal dotted line is the virial ratio of the singular 
isothermal sphere, (right) Virial ratio vs. radius calculated two different ways: (1) 
direct summation and (2) assuming spherical symmetry. 

halo results further by comparison with equilibrium halo models. 
2.4.1 Virial Ratio 

The scalar virial theorem states that for a self-gravitating system in static equilib- 
rium ((v) = 0) with no magnetic fields, 2T + W + S p = 0, where W is the potential 
energy, T is the thermal and kinetic energy, and S p is a surface pressure term, 

S p = - Jpv dS, (7) 

where dS is the surface area element. If the system is isolated, there can be no 
material outside to create a boundary pressure, and we have S p = 0, implying 
2T/\W\ = 1. Cosmological halos are not isolated systems, however, so we cannot 
expect S p = and 2T/\W\ = 1. In fact, the presence of infalling matter creates 
a flux of momentum into the halo which can act as a surface pressure (even for 
collisionless matter) in the virial theorem. With infall present, we therefore expect 
S P /\W\ < 0, implying 2T/\W\ > 1. 
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Figure 6: Virial ratio versus NFW concentration parameter for the halos as sim- 
ulated with and without gas, as they evolve. Solid, dashed, and long-dashed lines 
represent expected values for f3 = 0,0.5,1, respectively. Shown also are the values 
corresponding to the simulated halos, the isotropic singular isothermal sphere (SIS), 
and the truncated isothermal sphere (TIS). 



2.4.2 Singular Isothermal Sphere 

The singular isothermal sphere (SIS) is the singular solution of the equation of 
hydrostatic equilibrium with uniform temperature. The density is given by 

( \ = °o kBT (a\ 

P[r) 2^Gr 2 2^Gmr 2 ' { 1 

where do is the ID velocity dispersion, and T is the gas temperature. A more general 
class of solutions can be found, however, by allowing the anisotropy parameter (3 to 
be nonzero but remain independent of radius. The relevant equation is the Jeans 
equation [IB] . 

where the anisotropy parameter is now the "Eulerian" quantity, (5 = /3e> as defined 
in ^12.31 We can solve the Jeans equation in this case for the velocity dispersion 
as a function of (3 and the isotropic velocity dispersion ctq for the same mass 
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distribution as in equation (JHJ), to show 

2 3 — 2/3 2 

= sa^)" ' (10) 

If the anisotropy (3 does not vary with radius, we can determine the virial ratio 
of an SIS for a given value of the anisotropy, as follows. The potential energy at r is 



W(r) = J (pcb)dV = -2M(r)0fi. (11) 

The kinetic (or thermal) energy T is given by 

T=\J P (v 2 )dV = ^J po 2 dV = \M(r)a\ (12) 

Using equation (|l(Jj) to relate the actual a 2 to the one in the isotropic case a 2 , one 
obtains 

r = W^T) M < < 13 > 

The virial ratio is therefore 

2T = _3-2^ 

\W\ 2(1-/3) 1 ' 

2.4.3 Nonsingular Truncated Isothermal Sphere 

The Truncated Isothermal Sphere (TIS) model for the equilibrium halos which form 
from the collapse and virialization of top-hat density perturbations |128j (to be 
discussed here in fJQ) yields a virial ratio which differs significantly from the fa- 
miliar value of 2T/\W\ = 1 for isolated halos (i.e. with zero boundary pressure) 
with /3 = 0. The TIS model is the minimum energy solution of the isothermal 
Lane-Emden equation with non-singular, finite boundary condition for the central 
density. Since the TIS is a unique solution given by the minimum energy at fixed 
boundary pressure, the virial ratio is always the same value, namely 2T/\W\ — 1.37 
at (t = 29.4, where Ct is the truncation radius in units of the core radius. This 
value is smaller than that for the isotropic singular isothermal sphere (with nonzero 
boundary pressure), 2T/\W\ = 1.5, and is near the global minimum for all values 
of (, which is 2T/\W\ ~ 1.36 and occurs at ( ~ 22.6. At the intermediate radius 
£200 — 24.2, defined to be the radius within which the mean density is 200 times 
the background density, the TIS virial ratio has a value 2T/\W\ ~ 1.36. As seen in 
Fig. |3 the inner core region of the TIS is dominated by kinetic (or thermal, in the 
gas case) energy, whereas the value approaches that of the SIS at large (, where the 
core region becomes small compared with the size of the TIS and the density profile 
asymptotically approaches that for the SIS (p oc r -2 ). 



The Equilibrium Structure of CDM Halos 



14 



2.4.4 NFW Halos 

The equilibrium structure of halos with an NFW density profile was investigated by 
|92| . Using different values of /3(r), they found several analytical solutions to the 
velocity dispersion of the halo by integrating the Jeans equation for a given p(r) and 
P(r). In order to integrate the Jeans equation to find the velocity dispersion, it is 
necessary to set the velocity dispersion at some r. In the absence of a physical value 
for o r at the boundary of the halo, as could be inferred from some infall solution, 
the only other reasonable choice is to have u T — ► as r — > oo. The results obtained 
by fixing o r = at infinity should reflect the general trends associated with a single 
halo evolving in a quasi-equilibrium state. 

Shown in Fig. HO is the virial ratio versus concentration parameter for different 
values of 0(r) = (3q, as expected from the Jeans equation. The more isotropic the 
NFW halo, the lower the virial ratio. This is consistent with the fact that the 
surface pressure term is directly related to the radial velocity dispersion. The more 
concentrated the NFW halo, the lower the virial ratio. This is expected, since the 
more isolated a halo, the less relative importance of the boundary pressure, implying 
2T/\W\ -> 1. 

2.4.5 Simulated Pancake Instability Halos 

Virial ratios 2T/\W\ were calculated for both simulation runs, with and without 
gas. For the case with no gas, 

T(r) = (15) 

i 

where the sum is over all particles within r. For the simulation with gas included, 

T ( r ) = Y \ m i v 1 + Y 2 mikBTi - ( 16 ) 

i i 

The potential energy was found using the assumption of spherical symmetry, 

W{R) = -Y, r ^h, (IT) 

i n 

where Mi is the mass interior to rj and fi = /(rj) is a function which represents the 
particular form of the softening used in the P 3 M algorithm, and the sum is over all 
particles within R. As shown in Fig. El the assumption of spherical symmetry yields 
results which are very close at R = r2oo to those arrived at from the more rigorous 
definition, 

^) = -5EE- ( 18 ) 



2 ; 



where r {j = \n - rj| and fcj = /(r^). 
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Figure 7: Mass accretion history of simulation with dark matter only (dashed) 
along with the best-fit functional form from CDM simulations by |15L)j (solid) with 
5 = 2 and Of = 2.5a c (see equation Curve labeled "Gadget" was for the same 
initial conditions, run with the publicly available code of the same name |136| . 

Shown in Fig. are the virial ratios of both halos, with and without gas, plotted 
versus their concentration parameter. Each point corresponds to a different time in 
the halo's evolution, since the concentration of each halo changes with time. While 
the simulated halos lie below the expected curve for NFW halos of comparable 
anisotropy, the trend of smaller virial ratios for more concentrated halos is clearly 
evident, though there is significant scatter, particularly in the simulation without 
gas included. 

2.5 Halo Evolution 

2.5.1 Mass Accretion History 

The mass growth of the halo proceeds in three stages (Figs. and IHJ). Before 
a/a c ~ 3, the mass within r2oo grows very quickly, indicating initial collapse of the 
central overdensity. After a/a c ~ 3, the infall rate drops, and can be well-described 
by M oc a. Such an infall rate is reminiscent of self-similar spherical infall onto 
a point mass (e.g. JH). This infall rate cannot persist indefinitely, since there is 
only a finite mass supply to accrete onto the halo because of the periodic boundary 
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Figure 8: Evolution of the dark matter only halo in the intermediate (self-similar) 
regime, (left) Virial radius in units of expected caustic radius as defined in the 
text, and the average in the range 3 < a/a c < 7.5 (horizontal dotted line), (right) 
Halo mass and the linear best-fit (dotted). The vertical dotted lines indicate the 
virialization epoch at a = 3a c . 



conditions. As a consequence, the accretion rate slows after a/a c = 7. This is also 
expected to occur in halos forming from more realistic initial conditions, given that 
neighboring density peaks of a similar mass scale prevent any one halo from having 
the infinite mass supply necessary to sustain this mass accretion rate. In fact, as 
shown in Fig. the halo can be fit at nearly all times, especially later, by the more 
general fitting function 

M(a) = Moo exp [Sat/a] , (19) 

where S = [d In M/d In a] af is the logarithmic slope at a = af. For S = 2, we find 
a best-fit value of af = 2.5a c . This form was first used to fit the evolution of halos 
formed in CDM simulations f |15U| ^ We have identified three distinct phases in the 
halo evolution: initial collapse, steady infall, and infall truncation due to finite mass 
supply. In realistic collapse, we see that the halo evolves continuously from one 
stage to the next as evidenced by the continuous change in logarithmic slope of the 
fitting function, given by 

^ = S«. (20) 
a In a a 

In the intermediate stage of collapse, the resemblance to self-similar spherical 
infall onto a point-mass perturbation is evident in (Fig. |S1 and [^J) . Plotted in Fig- EI 
is the dark matter radial dimensionless velocity profile V(A) for different a/a c as 
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Figure 9: Radial velocity profiles in dimensionless units as in self-similar spherical 
collapse. The thick solid line is the radial velocity profile for an ideal 7 = 5/3 gas 
as in |14j . 



simulated with and without gas. The dimensionless velocity is 

t 2 / a \ 1 / 6 v r 

r ta {t) 3H 7 \a 7 J r taJ 

where we have used the relations r ta oc t 8 / 9 and a oc r ta is the turnaround 
radius, aj = 7a c , r ta ,7 = r ta (aj) is set by finding the radius at which v r = at 
a = aj, and the Hubble constant Hj = H(a = a 7 ). The dimensionless radius is 

r r / a \ ~ 4 / 3 



A = — = . (22) 

rta, r ta j \7a c J 

In the case of self-similar infall, this profile does not change with time. As seen in 
Figure IHJ V(A) for the simulated halo follows the self-similar solution closely, with 
A20o/A c — 0.78 (see Fig. |HJ, where A c is the radius at the outermost caustic and is 
approximately where the shock occurs in the collisional solution, and A200 = ^20o/^ta- 



2.5.2 Density Profile 

Although the halo generally grows by self-similar accretion for 3 < a/a c < 7, the 
mass density is better fit by an NFW profile with only one free parameter, c. Shown 
in Fig. E3 is the concentration parameter versus scale factor. We find it can be 
well-fitted by 

c = c f -, (23) 

Of 
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Figure 10: Evolution of concentration parameter of the best-fitting NFW profile for 
the dark matter halo in simulations with (right) and without (left) gas. Dotted line is 
the actual evolution, with the filled squares showing the mean concentration binned 
in scale factor. Errorbars indicate RMS fluctuations within each bin. The solid lines 
are the best-fitting linear evolution given by cnfw = Pf(a/af), with Of = 3a c . 

where af is the scale factor at which the accretion rate becomes proportional to a, 
marking the end of the collapse phase. The value df = 3a c is used here for both 
simulation cases and corresponds to the vertical dotted line in Fig. The solid lines 
in Fig. EI] correspond to the best-fit values Cf = 4.3 and Cf = 3.8 for the cases with 
and without gas included, respectively, where each data point was weighted by the 
goodness of the corresponding NFW profile fit. 

This linear evolution was also found by |22j and |15flj . In the latter case, they fol- 
lowed the mass accretion and merger histories of individual halos in a high-resolution 
CDM simulation of halos in the mass range ~ 10 11 — 1O 12 M . The mass accretion 
histories allowed them to determine a collapse epoch for each halo, which they cor- 
related with the halo's concentration. They found a best-fit slope for the linear 
evolution of concentration parameter of Cf = 4.1. We conclude that such an evo- 
lution of mass and concentration is generic and not limited to halos forming from 
CDM initial conditions. 
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3 The Fluid Approximation: ID Halo Formation From 
Cosmological Infall 

3.1 Self-Similar Gravitational Collapse 

3.1.1 Previous analytical models for halo formation: the spherical infall 
model 

Analytical approximations have been developed to model the formation of halos by 
the ID growth of spherical cosmological density perturbations, in either a collision- 
less gas or a fluid. We shall need to refer to some of these solutions to justify our 
fluid approximation in ^13. 21 It is necessary to begin, therefore, by briefly recounting 
this earlier work. 

[61] first presented the concept of the so called "secondary infall model (SIM)". 
This SIM refers to the effect of the addition of a point mass to a uniform, expanding 
Friedmann-Robertson- Walker universe as a perturbation which causes the surround- 
ing spherical shells to decelerate relative to the background universe, until they reach 
a radius of maximum expansion and recollapse. Subsequent work generalized this 
approach to include spherically-symmetric initial perturbations for which the over- 
density profile depends upon radius or mass as a scale-free power-law. Along this 
line, [IS] studied the dynamics of collisionless CDM halos using a self-similar model, 
adopting a scale-free initial overdensity parametrized by its shape: e in equation 
l|24j). [Tl] studied a special case of 0S] but also extended the analysis to a colli- 
sional fluid. showed that a power-law power spectrum would indeed generate a 
scale- free initial condition, such as was adopted by [IS]. They then argued that the 
resulting nonlinear structure would be described by a power-law profile determined 
by the shape of the power spectrum. 

Previously mentioned works adopted a rather unrealistic condition for the col- 
lisionless case, that of purely radial motion. N-body simulations of CDM halo 
formation find that the virialized region tends toward isotropic random velocities. 
Some attempts to incorporate tangential velocities within the framework of spher- 
ical symmetry have also been made by |9l 1651 1122"] . Along this line, one may also 
refer to work by [401 I63[ I89[ I91j and references therein. 

The fluid approximation for collisionless CDM halo formation has emerged re- 
cently. |142j showed that one could use fluid-like conservation equations to mimic 
the radial-only SIM models such as the [IS] model. [137] extended this analysis to 
include tangential motion. The fluid approximation has been used in the literature 
of stellar dynamics, but its application to CDM halo dynamics is rather new, and 
as will be described in H3.2[ it simplifies the description of dark matter dynamics 
substantially. 

Self-similar spherical infall models have also been used to study the effect of 
incorporating additional baryonic physics. P, ^H] and [110] have studied the effect 
of gas cooling on galaxy formation using self-similar models. Because the system 
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was tuned to maintain self-similarity in the presence of cooling, the cooling function 
used is not perfectly physical. However, as shown by P, one can use these models to 
test one's hydrodynamic code in the presence of cooling. And these models capture 
the generic behavior of galactic dynamics in the presence of realistic cooling. 

Our model [2113]) which will be described in §E3J is the first self-similar model to 
include the effective heat conduction resulting from SIDM collisionality. It utilizes 
the fluid approximation to describe SIDM halo dynamics, justifying its validity 
rigorously. The resulting formalism is similar to that of self-similar cooling models 
such as However, the focus is on dark matter physics rather than on baryonic 
physics. Thus, this is the first self-similar model in the presence of dark-matter-based 
heat conduction. It may also be used as a testbed for hydrodynamic simulations 
incorporating heat conduction. 



3.1.2 Halo formation from scale-free linear perturbations 

In the Einstein-de Sitter (EdS) background universe, an initial linear perturbation 
whose mass profile is spherically symmetric and has a scale-free, power-law form 

results in structure formation which is self-similar ( |48j ) . Each spherical mass shell 
around the center expands until it reaches a maximum radius (turnaround radius 
rta), and recollapses. For a given e, we have 

r ta oc t*, (25) 

where 

(48 J. Since there are no characteristic length or time scales for this problem other 
than the turn-around radius r ta and the Hubble time t, the gravitational collapse 
which ensues from this scale-free initial condition must be self-similar as long as the 
background universe is Einstein-de Sitter, in the absence of physical processes which 
introduce additional scales (e.g. SIDM collisionality). 

In general, if the unperturbed matter is a cold fluid, the infall which results 
from this perturbation is highly supersonic and is terminated by a strong accretion 
shock which thermalizes the kinetic energy of collapse. The accretion shock radius 
is guaranteed by self-similarity to be a fixed fraction of rt a (t) at all times. The 
mean density of the postshock region is, therefore, always a fixed multiple of the 
cosmic mean matter density. For most cases of interest, this postshock region is 
close to hydrostatic. For a collisionless gas, a similar description applies as long as 
the infalling matter initially had small (or no) random motions. In that case, each 
mass shell collapses supersonically as a single stream until it encounters a region of 
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Figure 11: Self-similar collisionless halo formation for e = 1: Comparison of the 
skewless- fluid approximation to the exact collisionless solution by ^3]. Solid lines 
represent the solution obtained from fluid approximation in the radial direction, 
while dotted lines represent the collisionless Bertschinger solution. Spikes in the 
density plot simply represent infinite values, corresponding to caustics, and therefore 
there is no physical significance in the height of these spikes. However, spikes in the 
velocity plot are finite and real. Note that solid lines do not represent the 7 = 5/3 
fluid Bertschinger solution. 



shell-crossing and density caustics, which encompasses all previously collapsed (i.e. 
interior) mass shells. All collapsed mass shells inside this region oscillate about the 
center. The radius of this region of shell-crossing, given by the outermost density 
caustic, is analogous to the shock radius in the fluid case. 

Results for the purely collisionless case were presented for several values of e by 
|48j and for e = 1 by ^3] (where the latter included a fluid component, as well). 
Figures ^2 and El show the exact similarity solutions we have derived ( |Hj ) for the 
purely collisionless cases with e = 1 and e = 1/6, respectively. As we describe below 
in § I3.1,3( these values roughly bracket the range relevant to cosmological halos in 
a CDM universe. We will refer to these solutions again for comparison in deriving 
our fluid approximation in § 13.21 
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Figure 12: Same as Fig. I1H but e = 1/6. Note again that solid lines were not 
generated from the 7 = 5/3 fluid approximation, but rather from the radial-only 
fluid approximation). 

3.1.3 Halo formation from peaks of the Gaussian random noise primor- 
dial density fluctuations 

The theory of halo formation from peaks in the density field which results from 
Gaussian-random-noise initial density fluctuations draws an interesting connection 
between the average density profile around these peaks and the shape of the fluctu- 
ation power spectrum. According to |64j . local maxima of Gaussian random fluctu- 
ations in the density can serve as the progenitors of cosmological structures. They 
show that high density peaks (y > 3, where v corresponds to vom peak) have a 
simple power-law profile 2 

A (r) cxr-( n+3 \ (27) 

where Ao(r) is the accumulated overdensity inside radius r, and n is the effective 
index of the power spectrum P(k) approximated as a power-law P(k) oc k n at 

2 |11| also get a similar result: local density maxima have a triaxial density profile, but as v 
increases it becomes more and more spherical with a profile converging to equation 1271 . 
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Figure 13: Effective index of the power spectrum (P(k) oc k n ) vs. halo mass for the 
ACDM universe. 

wavenumber k which corresponds to the halo mass through 



where po is the present cosmic mean matter density and R is the comoving, spatial 
top-hat filter radius 3 . The overdensity Ao(r) is equivalent to the fractional mass 
perturbation 5M/M inside radius r, 



From equations (|24|) and (|29j) . we deduce that the power-law power spectrum nat- 
urally generates a scale-free initial condition with 



According to this model, halos of a given mass M originate from density pertur- 
bations given by equation (|29|) with n determined by the primordial power spectrum 
after it is transferred according to the parameters of the background universe and 
the nature of the dark matter. We plot this effective n as a function of halo mass 
in Fig. EH for the currently-favored ACDM universe. The value of n ~ —2.5 is a 
reasonable approximation for galactic halos [i.e. n ~ —2.5 ±0.1 for M ~ 1O 8±2 M , 
while n ~ —2.5 ± 0.2 for M in the range from 10 3 M Q to lO n M ]. For halos in the 
cluster mass range, M ~ 1O 15 M0, n ~ —1.5. 

3 As long as one is concerned about the average profile, equation 12711 holds for any value of v. 
For small v, however, random dispersion around this average profile becomes substantial, limiting 
the generality of equation 1271 . 




(28) 



7T. + 3 

A (r) = 5M/M ocM~—. 



(29) 



e = (n + 3)/3. 



(30) 
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3.2 Fluid Approximation of Collisionless CDM Halo Dynamics 

This approach is used to simplify our understanding of the formation of CDM halos. 
We show that the fluid conservation equations for a gas with adiabatic index 7 = 5/3 
are a good approximation to the dynamics of CDM halos and self-interacting dark 
matter halos, which will be treated in detail in ^3.41 This may bother some readers 
since, strictly speaking, the collisionless nature of CDM prohibits the use of such an 
approximation. However, a couple of simple assumptions enable us to treat CDM 
halo dynamics with the usual fluid conservation equations. For a more detailed 
description, see 

We first assume spherical symmetry, and define the average of physical quantities 
as follows: 

P = I fd 3 v, (31) 



u = (v r ), (33) 

Pr = p{(Vr-(v r )) 2 }, (34) 

Po = P((v e - (v e )) 2 ) = p{vl), (35) 
V4> = p{(v<t, - {v<t>)f) = p(vl), (36) 

where / is the distribution function defined such that f(r, v)d 3 rd 3 i> = mass within 
an infinitesimal volume d 3 rd?v at (r, v), p is the density, (A) is the average value 
of a certain quantity A, u is the radial bulk velocity, p r is the "effective radial 
pressure", and pq is the "effective tangential pressure". Note that (vq) = {v^) = 
and pq = p^ because of spherical symmetry. Anisotropy in the velocity dispersion 
occurs in general for collisionless systems - i.e. p r ^ pe, or anisotropy parameter 
f3 ^ 0, where f3 = 1 — — implying that p r and pg should be treated separately. In 
a highly collisional system, which is well described by fluid conservation equations, 
p r = pg and the usual pressure p = p r = pg. 

A self-gravitating system of collisionless particles can be described by the colli- 
sionless Boltzmann equation 

! = g +v . V/ -V..g-0. (37) 
In spherical symmetry, / = /(|r|, v), equation (|37|) becomes 

= 9f | Vr 9f | ( V 8 +V l 

dt dr \ r 

-^(v r + v e cote)^-, (38) 




r 
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where $ satisfies the Poisson equation, V 2 <J> = AirGp |16j . By multiplying equation 
(|38|) by v™Vg, where m, n are integers, and integrating over velocity d?v, we can 
form a set of moment equations. Moment equations from the lowest order are 

| + ^ (r2(/5M)) = ' (39) 



d , . d . 2 , 2 o Gm 



<£&) + »£- r " (41> 



where m is the mass enclosed by a shell at radius r, = + ujp, and 

ri = £ (2 (tv - («r»«?) - (- 2 p(K - K» 3 )) , (43) 

^ = -^ r ^ P {(v r -{v r ))vl)). (44) 

Equations (|H9B42j) are conservation equations of mass, momentum, "radial" energy, 
and angular momentum, respectively. Note that equations (|39|) - (j44|) are all in exact 
form, and the hierarchy of equations is not closed in principle. 

Now we make a further simplification that the distribution of v r is skewless - 
vg and are naturally skewless because of spherical symmetry. In other words, 
we assume that v r has a symmetric distribution around (v r ). It is not straightfor- 
ward to show that Ti and T2 are negligible in equations (|41jl and (|42jl. However, we 
demonstrate that the assumption of "skewlessness" in the fluid approximation yields 
results which are in good agreement with the purely collisionless CDM structure for 
specific examples. Equations (jBHl) - (|4*T|) with the condition pg = and Ti = 
can be used to solve purely radial problems, such as the self-similar spherical infall 
problems with similarity solutions by Jl] (e = 1) and [48_. (e = 1/6). Comparisons 
of the fluid approximation in these shown in Figures ITTI and WR reveal an 

excellent agreement with the true collisionless solutions. The difference observed at 
caustics - places where the density becomes infinite - is negligible, because caustics 
do not affect the overall dynamics of the halo. Since the skew-free assumption nat- 
urally neglects dynamically unimportant structure (e.g. caustics) while accurately 
reproducing the profile of the exact solution in these radial cases, it may also be 
applied to describe CDM halos, in which particles also have a tangential motion. 

The final assumption is that the velocity distribution is isotropic, or p r = pg. 
This is an empirical assumption: CDM halos in cosmological N-body simulations 
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show mild anisotropy. For instance, p?j show that CDM halos in their numerical 
simulation can be well-fitted by a fitting formula 

P(r)=P m ^—, (45) 
r z + 4 

where r is in units of r2oo, and (5 m ~ 0.5 in N-body simulations |271 I31| . 

With these assumptions - spherical symmetry, skew-free velocity distribution 
and isotropic velocity dispersion - and with the operation (eauationl41|) + 2 x (equation 
HU, the usual energy conservation equations are obtained. They are 

| + ^ (r2(pn)) = °' (46) 
d , d , 9n 2 o Gm 



^{pu) + j^iv + pu 2 ) + -pu 2 = -p^f, I 17! 

777^1 ~ -Z-£r:(r*u), (48) 



D 3p p d 
D~V2~p> ~ ~~pr^dr' 

which are identical to the fluid conservation equations for a 7 = 5/3 gas in spherical 
symmetry. 

In this section, we showed how one could use the usual fluid conservation equa- 
tions to approximate collisionless CDM halo dynamics. ^13.31 and 13.41 describe its 
practical application. 



3.3 Halo Formation by Non-Self-Similar Infall: Mass Assembly 
History and the Origin of CDM N-body Halo Profiles 

N-body simulations of CDM have not only found a universal halo density profile, 
but have also found that the masses and concentrations of individual N-body CDM 
halos grow over time according to simple universal formulae |15L)j . Very similar 
results were also found for N-body simulations of halos formed by the instability of 
cosmological pancakes (see In what follows we use the fluid approximation to 
show that this universal time-dependent halo density profile can be understood as 
the dynamical outcome of continuous infall according to the universal mass accretion 
history [3]. 



3.3.1 Models and Initial Conditions 

We attempt to understand the form and evolution of dark matter halos with three 
spherically-symmetric models; halted infall, radial orbits, and a fluid approximation. 
Each model assumes that the mass M v i r within an overdensity A v i r follows the 
relation given by |15()j 

M vir (a) = Moo exp [Sa { /a] , (49) 
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Figure 14: (left) Density profile from halted infall model along with best-fitting 
NFW profile for this profile at present. Inset in upper-right shows same over much 
larger range, (right) Evolution of NFW concentration parameter in the halted infall 
model, compared with empirical relation of |15Uj for CDM N-body halos. Different 
line types indicate different ranges x- m < x < 1, within which halo was fit to an 
NFW profile, where x = r/r v i r , r v ; r = r2oo- 



where S is the logarithmic mass accretion rate dlriM v - n / dlna when a = cif . Here and 
in |15U| . 5 = 2. Such a relation is claimed to be a good fit to the evolution of halos 
of different masses and formation epochs. We use A v ; r = 200, so that the halo has a 
mass M200 and radius r*2oo- We have found an initial perturbation profile consistent 
with equation (|14j) (for EdS), 

"U^^ln^, (50) 
M M Moo 

where <5; depends on the initial scale factor ai, Of, and 5 v i r , and M is the unperturbed 
mass. The parameter b = 1 if pressure or shell crossing are not present outside of 
r v i r . If they are present outside the halo however, the initial perturbation is not 
guaranteed to lead to the correct mass accretion rate. In our radial orbit and 
fluid approximation calculations below, where shell crossing and pressure are indeed 
present outside of r v ; r , we have found that the resulting mass is close to that of 
equation (f4*9|) if b is allowed to vary as a fitting parameter (6 = 1(0.7) in the fluid 
approximation (radial orbits) calculations). 
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□ Radial Orbits 




1 10 



Figure 15: Evolution of mass for the radial orbits (top) and fluid approximation 
(bottom) simulations, compared with the empirical relation of |15(J| for CDM N- 
body results. Shown above each are the fractional deviations A = (M exac t — M)/M. 



3.3.2 Halted Infall Model 



In the simplest model, we have assumed that infalling shells come to an abrupt halt 
upon crossing into the halo, so that the velocity is zero for r < r v i r . The mass of 
the halo is 



4tt 



(51) 



where ~p is the cosmic mean mass density at that epoch. Mass continuity implies 
the density p vu just inside the virial radius is related to the rate of halo mass and 
radius increase according to 



~luT ~ 47rpvir7Vir ^o~- 



(52) 



Differentiating equation (|51|) and combining with equations (|49jl and (|52|). one ob- 
tains 



Pvir 



1 + 



3a 



i -i 



(53) 
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Figure 16: Density profile at the end of the radial orbit simulation. 



where po is the mean background density at a = 1. The virial radius is given by 

-Saf f 1 



?"vir,0 



aexp 



- - 1 

a 



(54) 



Equations (|53j) and (|54j) are parametric in a, implying a radial density profile p{r) = 
Pvir( r vir) that is frozen in place as matter crosses r v ; r . Taking the limit in which 
a — > oo, the outer density profile approaches p oc r -4 at late times, consistent with 
finite mass, while the inner slope becomes asymptotically flat. The NFW profile is 



p(x) ^ <5 vir ff(c) 
~p 3x(l + cx) 2 ' 



where 



9(c) 



(55) 
(56) 



ln(l + c) - c/(l + c) ' 

and x = r/r v i T . Combining equations 1)53(1 and (|55|) with x = 1, yields an equation 
for the evolution of concentration with scale factor (see Fig. I14|) . 



S* 



Of 



;i+c) 

9(c) 



(57) 
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3.3.3 Radial Orbits Model 

We use a finite-difference spherical mass shell code to follow the evolution of a small 
amplitude initial perturbation given by equation 1)50(1 . which is chosen so that the 
resulting virial mass will evolve according to equation (|49|). The shell code has an 
inner reflecting core and the results presented here used 20,000 shells. The resulting 
evolution of halo mass and the comparison of the halo density profile with the NFW 
profile are shown in Figures Hoi and [TBI 

3.3.4 Fluid Approximation Model 

As mentioned earlier in H3.21 the collisionless Boltzmann equation in spherical sym- 
metry yields fluid conservation equations (7 = 5/3) when random motions are 
isotropic. Halos in N-body simulations have somewhat radially-biased random mo- 
tion, but the bias is small, especially in the center. Outside the virialized halo, in 
the infall region, the radial bias is irrelevant, since the motion is highly supersonic 
and random motions do not affect the dynamics there. This isotropic fluid model is 
therefore a better approximation to halo formation in N-body simulations than one 
with purely radial motion. We use a 1-D, spherical, Lagrangian hydrodynamics code 
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Figure 18: Evolution of concentration parameter with scale factor in the isotropic 
fluid calculation. 

as in |143j . using 1,000 zones logarithmically spaced in mass. The initial conditions 
were chosen in the same way as those for the radial orbit model (Eq. (|50j0 . with 
zero initial temperature. Results are plotted in Figures 1151181 

3.3.5 Results 

Our results can be summarized as follows: 

• The halted infall model of ^3.3.21 does not reproduce the linear evolution of 
concentration parameter with scale factor reported by |150j . but can be fit by 
an NFW profile over a limited range of radii and scale factors (Fig. 1141) . 

• We have derived an initially linear perturbation profile that is a good model 
for the spherically-averaged initial condition that leads to the Lagrangian mass 
evolution of CDM halos found in N-body simulations by (150] . 

• Starting from this perturbation, the radial orbit model of section 13.3.31 fails to 
reproduce the inner slope of the NFW profile, approaching p oc r~ 2 instead, 
consistent with the argument of |120j (Fig. fIH)l . 
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• The fluid approximation of N3.3.41 however, leads to a halo that is well-fitted 
by the NFW and Moore profiles for all radii resolved by N-body simulations 
(r/r 200 > 0.01) (Fig. EJ). 

• In addition, the evolution of the NFW concentration parameter in the fluid 
approximation is a close match to that of |150j . with cnfw = 4.25a/af a good 
fit (Fig. H8|) : cnfw = 4.1a/cif was the relation reported by |150j . 

The fluid approximation model reproduces the N-body results remarkably well, 
once the mass accretion history is given. We are thus led to conclude that compli- 
cated merging processes are not necessary in order to understand the overall struc- 
ture and evolution of the halo mass distribution, and that it is largely determined 
by the mass accretion history. 

3.4 Structure of Self-Interacting Dark Matter Halos 

As we will demonstrate below, halo models with a small flat core, such as the 
Truncated Isothermal Sphere (TIS) halo equilibrium model described in § |1] agree 
well with the statistical properties of observed halos and with the averaged properties 
of simulated halos. However, the dynamical origin of such an equilibrium structure 
and the particular question of how the flat density core can form are not addressed by 
an equilibrium analysis. Halos in high-resolution N-body simulations in the ACDM 
model, on the other hand, show cuspy density profiles. 

Among several suggestions made to resolve this problem, self-interacting dark 
matter (SIDM) has drawn substantial attention. It was suggested by j!35j that the 
purely collisionless nature of CDM be replaced by SIDM which interacts by non- 
gravitational, microscopic interaction (e.g. elastic scattering). Various attempts 
have been made to incorporate SIDM in semi-analytical studies (e.g. []"""]) as well 
as in N-body simulations (e.g. [2*51 13*51 1551 1 154j ). and they all show that flat density 
cores can arise. However, new theoretical problems have arisen for these SIDM 
halos: 1) study of "isolated" (i.e. no cosmological infall) halos revealed unstable 
cores ( |25[ I85|) and 2) cosmological N-body simulations ( |36[ 1154] ) showed a more 
prominent flattening effect for high-mass halos (M halo ~ 10 12 M Q ) than for low -mass 
halos (M halo ~ 1O 9 -1O 1O M ), which appears to contradict observations which show 
that the flattening effect is most prominent for dwarfs and low surface brightness 
(LSB) galaxies with M halo ~ 10 9 - 1O 1O M . 

The semi-analytical study by 10 applied the fluid approximation to an isolated 
halo, properly accounting for the effective "conduction" arising from the SIDM 
collisionality. They showed that SIDM halos could have a lifetime long enough to 
survive in a Hubble time, if initially the scattering mean free path is much greater 
than the size of the halos. However, as was found in |25| I85j . the scattering mean 
free path can be smaller than the halo size if the scattering cross-section is large 
enough. In this case, isolated SIDM halos would still be unstable. 
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Figure 19: Halo density profile produced by fluid approximation with e = 1/6 for 
standard CDM halos, compared to the best-fitting NFW and Moore profiles. The 
adiabatic solution, whose inner slope is about -1.27, is a good fit to density profiles 
of CDM halos from N-body simulations. 

We show that these problems can be naturally resolved [21 El- We partially 
adopt the formalism by ^U] and improve upon it by properly including cosmological 
infall. We have thereby derived an analytical, cosmological similarity solution. Our 
approach is as follows. As shown in § 13.1.31 the power spectrum in the mass range 
of dwarfs and LSBs is well described by a power law P(k) oc k n with n ~ —2.5. 
From equation ((2*9")) . we then have 

5M/M oc M -1 / 6 . (58) 

This also results in 

rta oc i 2 , (59) 

where r ta refers to a "turn-around" radius at which radial velocity vanishes. For 
any value of n, r ta serves as a natural length scale for self-similar accretion. 

This is not true for SIDM halos because the SIDM interaction brings in a new 
length scale which does not in general preserve this self-similarity. However, n = 
—2.5 is a magic number which preserves the self-similarity because in this case the 
new length scale is always a fixed fraction of the halo size. Specifically, SIDM 
introduces a conductive heat flux ( 10 ) 
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Figure 20: Profiles of the low-Q solutions. 

where a is the interaction cross-section, p is the density, and p = p (v — (v)) 2 is the 
effective "pressure" (HHH1 considered an elastic scattering case, in which a = 2.26 
and b = 1.002). The thermal energy changes according to oc r^ a t~ 5 , while the 
conductive heating rate is V • f oc rj^ a t~ 7 . Self-similarity is preserved only when these 
terms have the same time dependence, or rt a oc t 2 . From eq (|25|). this condition is 
equivalent to n = —2.5. 

As shown in ^3.1.31 dwarfs and LSBs are well described by self-similar solutions 
with this requirement. This implies that the core of an SIDM halo in that mass 
range would be stable, since the core's size would grow in proportion to the virial 
radius r v , or the turnaround radius rt a - In other words, the collapse of halo cores 
predicted for isolated SIDM halos is entirely prevented by cosmological infall, which 
resolves the issue of stability of SIDM halo cores. 

Before treating soft-core solutions, we describe properties of this self-similar 
solution without SIDM collisionality, i.e. the case for / = 0. The solution obtained 
by applying the fluid approximation - equations Q46I - H5|) - to the e = 1/6 (n = 
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Figure 21: Profiles of the high-Q solutions. 



—2.5) case yields a profile which resembles CDM halos in many aspects. First, 
as shown in Figure EE it has a density cusp with a logarithmic slope ~ —1.27 at 
4 x 1CT 3 < r/V 2 oo < 1-4 x 1CT 2 , which is between -1 (NFW profile) and —1.5 
(Moore profile) 4 . Second, the temperature profile is very similar to that of CDM 
halos. The temperature is zero at the center, rises to a maximum as radius increases 
until a point where the temperature starts to drop again (see the solid line in the 
temperature plot of Figure l2fl|) . This validates the use of our model for the study of 
SIDM halos and even makes it possible to correctly calculate the physical effect of 
SIDM collisionality. 

In the presence of SIDM collisionality, similarity solutions reveal soft-cores. Dif- 
ferent solutions arise for different values of the dimensionless collisionality parameter 
Q = op\,r s , where pb is the cosmic mean matter density and r s is the effective shock 

4 This solution is well-fitted by NFW profiles with very low concentration parameter c ~ 3. 
e = 1/6 corresonds to a fast mass accretion rate, din M / din a = 6, which is typical for the earliest 
epoch of mass accretion for standard CDM halos at a/a c < 1. Therefore, this solution can be said 
to capture the very early formation epoch of CDM halos. 
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Figure 22: Collisionality parameter Q vs. halo mass M. (left) "Typical" halos, 
well represented statistically by a 1 — om fluctuation. Different curves correspond 
to different scattering cross section a {a = 0.56, 5.6, 200, 1.2 x 10 4 , 2.7 x 10 4 from 
bottom to top), (right) Different v — cm fluctuations for a = 200cm 2 g^ 1 . 

radius. Q is constant during the matter-dominated epoch, z > 1, because pb oc t~ 2 
and r v oct 2 . 

Fig-EHJand Fig.[^show profiles for different Q's. One would naively expect that 
as Q (a) increases, the central density decreases because of the increasing flattening 
effect. However, if Q is too large, the mean free path between collisions is smaller 
than the size of the halo, which reduces the flattening effect. Therefore, there are 
two opposite regimes, low-Q (long mean free path; Fig |20j) and high-Q (short mean 
free path; Fig 121)) . which are separated by some threshold Q t h- We found that 
Qth = 7.35 x 10" 4 . 

The behavior of high-Q solutions also gives us a clue as to how a constant cross 
section a may still resolve the problem 2). By using the Press-Schechter formalism, 
E] relate Q to the halo mass M. When one matches Q(M ~ 10 9 - 10 10 M o ) 
to Qth) which is well-fitted by the empirical rotation fitting formula for dwarfs and 
LSBs (described in the following paragraph), one finds that a ~ 200cm 2 g~ 1 and 
collisionality increases as mass increases, which indicates that halos more massive 
than dwarfs and LSBs will have less density-flattening. This argument is consistent 
with previous N-body simulations by |154j and they used a ~ 0.1 — 10cm 2 g _1 , 
which is in the low-Q regime. In the low Q regime, increasing Q means more density- 
flattening. They simply did not push a to higher values, which would have covered 
the high-Q regime (See Fig. [2*2*)l . 

Since the Qth solution shows the most prominent density-flattening effect, it can 
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Figure 23: Rotation curve fitting: The Qth solution is the one that is best-fitted by 
the Burkert profile. 

be said to be the most relaxed state possible for SIDM halos. We found that in fact 
the solution with Q = Q t ^ is well described by the TIS profile [3]. We fitted rotation 
curves of various Q solutions with the empirical rotation curve fitting formula found 
by Burkert (|24|: Burkert profile) and found that the Q = profile is best fit 
by the Burkert profile (See Fig. I23|) . As shown in [JJ, the TIS profile is also well- 
fitted by the Burkert profile. The fact that Qth is best fit by the Burkert profile 
(and therefore the TIS profile) suggests that in fact Qth is the most relaxed system 
among the various Q solutions, since the requirement for TIS, the minimum energy, 
is usually met by the most relaxed system. 

4 The Truncated Isothermal Sphere (TIS) Model 

We have developed an analytical model for the postcollapse equilibrium structure 
of virialized objects which condense out of a cosmological background universe, 
either matter-dominated or flat with a cosmological constant [1281 171j . The model 
is based upon the assumption that cosmological halos form from the collapse and 
virialization of "top-hat" density perturbations and are spherical, isotropic, and 
isothermal. This leads to a unique, nonsingular TIS, the minimum-energy solution 
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Table 1: A comparison of three approximations for the postcollapse equilibrium 
structure of top-hat density perturbations in a low density universe. 



SUS a SlS b TIS (SCDM;ACDM, z = 0) 



v/vsvs 


1 


0.833 


1.11;1.07 


r/Tsus 


1 


3 


2.16;2.19 


Po/pt 


1 


oo 


514;530 


(p)/pt 


1 


3 


3.73;3.68 


n/r 


- NA- 


oo 


29.4;30.04 


A c /A CjSU s 


1 


1.728 


0.735;0.774 


K/\W\ 


0.5 


0.75 


0.683;0.690 



"A top-hat perturbation of a given mass collapses at a given redshift in a background universe 
with given values of fio and Ao; all of these values are held fixed in this comparison of the three 
approximations. 

These SIS numbers are an approximation which ignores the small modification of the Lane- 
Emden equation solution to take account of A 7^ 0, but accounts for the more important effects of 
A / on top-hat evolution, energy conservation and the virial theorem. 



of the Lane-Emden equation (suitably modified for non-zero cosmological constant 
A / 0). The size rt and velocity dispersion ay are unique functions of the mass 
and redshift of formation of the object for a given background universe. Our TIS 
density profile flattens to a constant central value, po, which is roughly proportional 
to the critical density of the universe at the epoch of collapse, with a small core 
radius tq r t /30 (where a v = AitGporQ and tq = rxing/3, for the "King radius" 
''King, defined by [TSj, p. 228). The density profiles for gas and dark matter are 
assumed to be the same (no bias), with gas temperature T = pm p a v /ks- While 
this TIS density profile is obtained numerically by solving a differential equation, it 
is well-fitted by the following approximation: 

p( r ) = Po ( 2 — 2T~2 ~ xa 2 ) > ( 61 ) 

\a +r z /rQ tr + r^/r^J 

where A = 21.38, B = 19.81, a = 3.01, b = 3.82 [T281I7T]. 

These TIS results differ from those of the more familiar approximations in which 
the virialized sphere resulting from a top-hat perturbation is assumed to be either 
the standard uniform sphere (SUS) or else a singular isothermal sphere (SIS). We 
summarize their comparison in Fig. [23] and Table 1, where rj is the final radius of 
the virialized sphere in units of the top-hat radius r m at maximum expansion (i.e. 
VSVS = 0.5), pt = p(rt), (p) is the average density of the virialized spheres, A c = 
(p)/Pcrit(^coii)) an d K/\W\ is the ratio of total kinetic (i.e. thermal) to gravitational 
potential energy of the spheres. 
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Figure 24: (top) Density profile of TIS which forms from the virialization of a top- 
hat density perturbation in a matter-dominated universe. Radius r is in units of r m 
- the top-hat radius at maximum expansion. Density p is in terms of the density 
Psus of the SUS approximation for the virialized, post-collapse top-hat. (bottom) 
Logarithmic slope of density profile. 



For all cases of current astronomical interest, an excellent approximation to the 
exact results in [7J (Paper II) for the dependence of the TIS halo parameters on the 
halo mass, collapse redshift and the background universe is given by [7J according 
to 

1 /3 

rm = 337 - 7 {wm^) [^oior 1 ^ 273 k P c, (62) 

1 /3 

n = 187 - 2 (l0^o) ^(^oiir 1 ^ 273 kpc, (63) 
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ro 



6.367 



M 



10 



12 



7.843 x 10 l 



M 
/i 



1/3 



[i^cou)] -1 /* -273 kpc 
M \ 2/3 



4 



1.098 x 10 



0.59 

M 



1O 12 M 

2/3 



1.799 
3.382 



1O 12 M _ 
10 4 [F(z coll )] 3 /j crit (z 



F(z con )h 2 / 3 K, 

F(z con )h 2 / 3 km 2 s- 2 , 
0) 



10- 25 [F(z coll )] 3 /i 2 gem" 3 . 



where 



^(^coll) = 



h{z co \\) 
h 



A 



c.TIS 



(^coll, Ao) 



A 



c,TIS 



(A = 0) 



n n a 



c,SUS 



Q(z co ll) 18^ 2 



1/3 



(1 + Zcoll) 



(64) 
(65) 
(66) 

(67) 
(68) 



For the EdS case, F = (1 + z co \\), while for an open, matter-dominated universe 

1 /3 

and a flat universe with a cosmological constant, F — > Q (1 + z C oii) at early times 
[i.e. x — > 0]. Here p is the mean molecular weight, where p = 0.59(1.22) for an 
ionized (neutral) gas of H and He with [iJe]/[If] = 0.08 by number. A C) sus is well- 
approximated by A Ci sus = 187r 2 +cix— C2X 2 , where x = fl(z co \\) — l, 0(z co u) = Oo(l+ 
zf [h/h(z)} 2 , where \h(z)/h} 2 = ft (l + ^) 3 + A (or Q (l + zf + (l-n )(l + z) 2 ) and 
ci = 82 (60) and ci = 39 (32) for the flat (open) cases, f^o + \) = 1 (^0 < 1> Ao = 0), 
respectively [2*T] . 



4.1 TIS Model vs. Numerical CDM Simulations 

The TIS model reproduces many of the average properties of the halos in numerical 
CDM simulations quite well, suggesting that it is a useful approximation for the 
halos which result from more realistic initial conditions: 

(1) The TIS mass profile agrees well with the fit to N-body simulations by |l()7j 
("NFW") (i.e. fractional deviation of ~ 20% or less) at all radii outside of a few 
TIS core radii (i.e. outside King radius or so), for NFW concentration parameters 
4 < cnfw < 7 (Fig. |23|) . The flat density core of the TIS halo differs from the 
singular cusp of the NFW profile at small radii, but this involves only a small 
fraction of the halo mass, thus not affecting their good agreement outside the core. 
As a result, the TIS central density po can be used to characterize the core density 
of cosmological halos, even if the latter have singular profiles like that of NFW, as 
long as we interpret po, in that case, as an average over the innermost region. 

(2) The TIS halo model predicts the internal structure of X-ray clusters found 
by gas-dynamical/N-body simulations of cluster formation in the CDM model. 
Our TIS model predictions, for example, agree astonishingly well with the mass- 
temperature and radius-temperature virial relations and integrated mass profiles 
derived empirically from the simulations of cluster formation by |46| I96j (EMN). 
Apparently, these simulation results are not sensitive to the discrepancy between 
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Figure 25: Profiles of density (left) and integrated mass (right), for TIS (solid) and 
NFW with cnfw = 4 (short-dashed), 5 (dotted) and 7 (long-dashed) with same 
(r 2 oo,M 2 oo)- 



our prediction of a small, finite density core and the N-body predictions of a den- 
sity cusp for clusters in CDM. Let X be the average overdensity inside radius r 
(in units of the cosmic mean density) X = (p(r))/p. The radius-temperature virial 

1 /2 

relation is defined as rx = rio(-X")(T/10keV) ' Mpc, and the mass-temperature 
virial relation by M x = Mi (X)(T/10keV) 1/2 h^W 15 M Q . A comparison between 
our predictions of rio(X) and the results of EMN is given in Fig. 1261 EMN ob- 
tain Mi (500) = 1.11 ± 0.16 and Mi (200) = 1.45, while our TIS solution yields 
Mi (500) = 1.11 and Mi (200) = 1.55. 

(3) The TIS halo model also successfully reproduces the mass - velocity dis- 
persion relation for clusters in CDM N-body simulations and its dependence on 
redshift for different background cosmologies. N-body simulation of the Hubble 
volume [(1000 Mpc) 3 ] by the Virgo Consortium jl^j yields the following empirical 
relation: 

r 1 1- 1 0.33 

av = (1080 ±65) h(z)M 2 oo/W 15 M Q km/s, (69) 

where M200 is the mass within a sphere with average density 200 times the cosmic 
mean density, and h(z) = foo\A^o(l + z)' 6 + ^0 is the redshift-dependent Hubble 
constant (assuming a flat background universe). Our TIS model predicts: 

oy = 1103 /i(z)M 2O o/lO 15 M km/s, (70) 

in excellent agreement with simulations. 

(4) The TIS model successfully predicts the average virial ratio, iif/|W|, of halos 
in CDM simulations. An equivalent TIS quantity, GM2oo/( r 200 cr y) = 2.176, is 
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Figure 26: Cluster radius-temperature virial relation for CDM simulation results 
(at z = 0) as fit by EMN (triangles) and as predicted by TIS (solid curve). 



plotted for dwarf galaxy minihalos at z = 9 in Fig. 1271 for simulations described in 
|125j and [22]) showing good agreement between TIS and N-body halo results. A 
similar plot, but of K/\W\ for such halos, was shown by [77j based upon N-body 
simulations, in which the average is close to 0.7, as predicted by the TIS 

model (Table 1). Those authors were apparently unaware of this TIS prediction since 
they compared their results with the SUS value of i^/|W|, 0.5, and interpreted the 
discrepancy incorrectly as an indication that their halos were not in equilibrium. 



4.2 TIS Model vs. Observed Halos 



(1) The TIS profile matches the mass profiles of dark-matter-dominated dwarf galax- 
ies deduced from their observed rotation curves . Equation ^ can be integrated 
to yield an analytical fitting formula for the TIS rotation curve, as well, given by 



v r 



cry 



A-B + 



bB tan 



-i 



-aAtan 



1/2 



(71) 



where £ = r/ro and v(r) = (GM[< r]/r) l l 2 . The observed rotation curves of dwarf 
galaxies are well fit by a density profile with a finite density core given by 



p(r) 



P0,B 



(r/r c + l)(r 2 /V2 + l) 
for which the rotation curve is given by 

v B (r) _ Jln[(CB + l) 2 (C| + l)] -2tan- 1 (Cs) 



1/2 



(72) 



(73) 
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Figure 27: Virial relation GM200 / ((Tyr2oo) predicted by the TIS model (horizontal 
line) vs. mass for halos from N-body simulations 125 [with lcr error bars]. 

where (b = f"A"o,B and u^b = (vrGpo,B^o b) 1 ^ 2 - ^he ^IS model gives a nearly perfect 
fit to this profile (Fig.l28|). with best fit parameters Po,b/pq,TIS = 1-216, rc/r^xis = 3.134. 
This best-fit TIS profile correctly predicts v m3 , x , the maximum rotation velocity, and 
the radius, r max , at which it occurs in the Burkert profile: r maX! B/Y maX) Tis = 1.13, 
and v max,B/^max,Tis = 1.01. (i.e. excellent agreement). 

(2) The TIS halo model can explain the mass profile with a flat density core mea- 
sured by |146j for cluster CL 0024+1654 at z = 0.39, using the strong gravitational 
lensing of background galaxies by the cluster to infer the cluster mass distribution 
|126j (see, however |34| for an alternative view on the structure of this cluster). The 
TIS model not only provides a good fit to the projected surface mass density distri- 
bution of this cluster within the arcs (Fig. I29|) . but also predicts the overall mass, 
and a cluster velocity dispersion in close agreement with the value o v = 1150 km/s 
measured by [IT] . 

4.3 Making Tracks on the Cosmic Virial Plane 

The TIS model yields (po> ^y,r t , ro) uniquely as functions of (M, z co \i). This defines 
a "cosmic virial plane" in (po, ro, ov)-space and determines halo size, mass, and 
collapse redshift for each point on the plane. In hierarchical clustering models like 
CDM, M is statistically correlated with z co \\. This determines the distribution of 
points on the cosmic virial plane. We can combine the TIS model with the Press- 
Schechter (PS) approximation for z co u(M) - typical collapse epoch for halo of mass 
M - to predict correlations of observed halo properties. 

According to the PS approximation, the fraction of matter in the universe 
which is condensed into objects of mass > M at a given epoch is / co ii(> M) = 



The Equilibrium Structure of CDM Halos 



44 



0.05 F 




Figure 28: Rotation Curve Fit. Solid line = Best fit TIS; Dashed line = Burkert 
profile. 

erfc 

(zv/2 1 / 2 ), where v = <5 crit /a(M), <r(M) is the standard deviation of the density 
fluctuations at that epoch, according to linear theory, when filtered on mass scale 
M, and 5 cr ;t is the fractional overdensity of a top-hat perturbation when this linear 
theory is extrapolated to the time of infinite collapse in the exact nonlinear solution. 
The "typical" collapse epoch for a given mass is that for which a{M) = 5 cr j t (i.e. 
v = 1). For a given z co u, this defines a typical mass scale: M* = M{y = 1). 

If we approximate the power-spectrum of density fluctuations at high redshift 
(e.g. just after recombination) as a power-law in wavenumber k, P(k) oc k n , and 
define a mass M oc A; -3 , then F oc M~ ( 3 + n )/ 6 if n = n e g = —3(2yp + 1), where 
UF = (cilni ? /(flnM)exact at the relevant mass scale. For all masses in the EdS 
case and for masses which collapse early in the low-density, matter-dominated and 
flat, universes, yp reduces to y a = (dlncr/dlnM) exact , and (1 + z co \\) oc a{M) oc 
M -(3+n)/6 ) ^ere a(M) is evaluated at the same cosmic time for all masses. The 
dependence of n e g and z co \\ on M for 1-a fluctuations is shown in Figure for 
ACDM, along with the approximate n e g which results if yp is replaced by y a , which 
shows that the latter is a very good approximation for all masses M < 10 12 M Q /i _1 . 

(1) The combined (TIS+PS) model explains the observed correlation of the 
maximum circular velocity v max and its location r max in the rotation curves of dwarf 
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Figure 29: Projected surface density of cluster CL 0024+1654 inferred from lensing 
measurements, together with the best-fit TIS model. 



spiral and LSB galaxies, with preference for the currently favored ACDM model with 
no tilt of the power spectrum of the primordial density fluctuations [Slj (Fig. 130ft . 
According to 103], the observed f max — ^max correlation can be expressed as follows: 

<W = ^(w/lkpcj^knis- 1 . (74) 

Our results for the ACDM case indicate that the galaxies which make up the v max — 
r max data points in Figure OH collapsed at redshifts 1 < z co \\ < 6 with masses in 
the range 8 x 10 9 < Mo/(M /i _1 ) < 3 x 10 11 . Hence, the precollapse fluctuation 
growth rate is approximately EdS, and we can let £l(z co \i) = 1- in that case, (1 + 
z co il) oc ct(M) oc M-( 3 + n V 6 , r max oc M^ +n V & % 1/3 and v max oc M( 1-n )/ :L2 fijj/ 6 , 
which combine to yield 

v — v (r Ir \( 1 -™)/[2( 5 + n )] 

"max — ''max,* \l max/ ' max,* J > \' °J 

where v m ax,* and r ma x,* ar e for a l-cr fluctuation of fiducial mass M*, with (l+2 C oil) 
' 1 + z Tec )a(M*, z rcc )/5 CT it, where <5 cr it = 1.6865 and <r(M*, z Tec ) is the value of cr(M*) 
evaluated at the epoch of recombination z rec [i.e. early enough that (1 + z)u is 
independent of z]. Over the relevant mass range M = lO lo±1 /i _1 M0, n e g ~ —2.4 ± 
0.1 for our ACDM case. For = 10 10 h~ 1 M Q , our COBE-normalized, flat ACDM 
case (Qq = 0-3, h = 0.65) yields (1 + z rec )a(M*, z rec ) = 5.563, so (1 + z C oil) = 3.30, 
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Figure 30: v max — r max correlation predicted by (TIS+PS) model for ACDM [COBE 
normalized, 0,q = 1 — Ao = 0.3, h = 0.65; no tilt, i.e. primordial power-spectrum 
index n p = 1), for halos formed from v — a fluctuations, as labelled with v- values. 
All curves are (TIS+PS) results, except curve labelled "Burkert" is a fit to data 
[2*1| . Observed dwarf galaxies (triangles) and LSB galaxies (squares) from [55] . 

"max.* = 53.2 km s ^ and r max * = 5.59 h x kpc. With these values and n e ff = —2.4, 
equation (f75]> yields the TIS model analytical prediction 

<W = (13.0 km s^Xw/l kpc) ' 65 , (76) 

remarkably close to the observed relation in equation (|74|) . 

(2) The TIS+PS model also predicts the correlations of central mass and phase- 
space densities, po and Q = po/ay, of dark matter halos with their velocity disper- 
sions ay, with data for low-redshift dwarf spheroidals to X-ray clusters again most 
consistent with ACDM with no tilt [TJZl (Fig.E2J). 5 A fully analytical approxima- 
tion for these correlations was also derived in |127j by combining the TIS+PS model 
with the power-law model for P{k) described above, which yields po oc M~( n+3 ^ 2 , 

5 There have also been claims that po =const for all cosmological halos, independent of their 
mass, as expected for certain types of SIDM |49l I78| . This claim, however, is not supported by 
most current data (Tig. 1321 . 
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Figure 31: (a)(upper) Effective logarithmic slope, n e s, of the density fluctuation 
power spectrum versus halo mass M = M+ (i.e. for l-cr fluctuations) for ACDM, 
based upon yp (solid curve) and the approximation which uses y a , instead (dashed), 
(b) (lower) Typical collapse redshift z co n for halos of mass M* in ACDM. 



oy oc M^ 1 "™)/ 12 , and Q oc AiH n+7 )/ 4 . These combine to give 

Q = Q*{w/crv,*) a ~ 3 , (77) 

and 

Po = p ,*(a v /(Tv,*) a , (78) 

where a = 6(n + 3)/(n — 1), except that a = for M > M+(z = 0) (e.g. for untilted 
ACDM M±(z = 0) w 1O 13 M ), for which z coVl = is assumed. 

Similarly, we can compare different flavors of the ACDM model, with no tilt 
(n p = 1), slight tilts (n p = 0.9 or 1.1) and with the "running" power-spectrum index 
suggested by the first year Wilkinson Microwave Anisotropy Probe results [134] 

n p (k) = 0.93 - 0.015 In ( OO5 M pc -0 ( Fi §- E3>- In this case the untilted ACDM 
model is again moderately preferred. 
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Figure 32: (a) (left) Maximum phase-space density Q = po/a v versus velocity 
dispersion ay for halos observed today, as predicted for various CDM universes by 
the TIS model combined with the Press-Schechter approximation (TIS + PS) (solid 
curves) for halos formed from v-a fluctuations, as labelled with the values of is, for 
v = 0.5, 1, 2, and 3. In each panel, lines representing halos of different mass which 
collapse at the same redshift are shown for the case z co \\ = 0, as labelled. Each panel 
represents different assumptions for the background universe and primordial den- 
sity fluctuations, as labelled: COBE-normalized ACDM (Ao = 0.7, f^o = 0.3) (up- 
per left), cluster-normalized SCDM (Qq = 1) (upper right), and COBE-normalized 
OCDM (Oo = 0.3) (lower panels), all assuming h = 0.7 and primordial power spec- 
trum index n p =1 (i.e. untilted), except for OCDM1.3, for which n p =1.3. Data 
points represent observed galaxies and clusters, taken from the following sources: 
(1) 49 late-type spirals of type Sc-Im and 7 dSph galaxies from (SUEZ! (open trian- 
gles); (2) Local Group dSph Leo I from [55] (filled square); (3) 28 nearby clusters, 
ay from [HH] and [20]) and po from [HZ] (crosses), (b) (right) Same as (a), except for 
po vs. ay. 
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Figure 33: Same as Fig.E21 but for ACDM variants: (top left) untilted ACDM (n p = 
1), (top right) ACDM with running power-spectrum index n p = n p {k) according to 
the best WMAP fit, (bottom left) ACDM with n p = 0.9, and (bottom right) ACDM 
with n„ = 1.1. 
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5 Gravitational lensing by CDM halos: Singular vs. 
nonsingular profiles 

5.1 INTRODUCTION 

The gravitational lensing of distant sources has in recent years become one of the 
most powerful tools in observational cosmology |13l 1135] . Since the effects of gravita- 
tional lensing depend upon the redshift of the source, the cosmological background, 
and the distribution of matter in the universe, they can be used to constrain the 
cosmological parameters and the primordial power spectrum of density fluctuations 
from which structure originates. In addition, many of the effects produced by grav- 
itational lenses, such as image multiplicity, separations, and time delay, depend 
strongly upon the matter distribution inside the lenses. Hence, measurements of 
these effects can provide a unique tool for probing the matter distribution inside 
collapsed objects like galaxies and clusters, providing the only direct measurement 
of their dark matter content, and constraining the theory of their formation and 
evolution. 

Until recently, the internal structure of halos adopted in lensing studies was 
generally some gravitational equilibrium distribution, either singular or nonsin- 
gular (e.g., King model, singular isothermal sphere, pseudo- isothermal sphere), 
not necessarily motivated directly by the theory of cosmological halo formation 

[m Ea cza E3 lEa nm im una esi eh ma unsj. a s the theory of haio for- 
mation in the CDM model has advanced in recent years, however, the halo mass 
profiles adopted for lensing models have been refined to reflect this theory. Nu- 
merical simulations of large-scale structure formation in Cold Dark Matter (CDM) 
universes predict that galaxies and clusters have a singular density profile which 
approaches a power law p cx r~ n at the center, with the exponent n ranging from 1 
to 1.5 r5lH f54| f55l r58^ [67^ rTTTOl rTT)2l rTTT7l rTTTHl f8T] rTT2l rT^4] . These results are 
in apparent conflict with observations of rotation curves of dark-matter-dominated 
dwarf galaxies and low surface brightness galaxies, which favor a flat-density core 
[2l)ll§!IHTUUllll8| . On the scale of clusters of galaxies, observations of strong gravita- 
tional lensing of background galaxies by foreground clusters also favor the presence 
of a finite-density core in the centers of clusters |146| l5T| 1123] . 

Lensing by the two kinds of halo mass profiles, singular versus flat-core, will be 
different. This has led to attempts to predict the differences expected if the halos 
have the singular cusp of the Navarro- Frenk- White (NFW) or Moore profiles or else 
a profile with a flat core [ZH1 EH1 E01 EZEJ EDI EEH1- Several authors have studied 
the effect of lensing by halos with a flat-density core [H E3 ESI EEH HHD EES] or 
by NFW or Moore profiles that have been generalized, so that the inner slope of 
the density profile is arbitrary [7^1 IHB 11211 1153] . These particular density profiles 
are essentially mathematical conveniences without physical motivation. There is no 
underlying theoretical model in these cases that was used to predict the value of the 
core radius or the departure of the inner slope of the density profile from the value 
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found by N-body simulations of CDM. By contrast, the TIS model is based on a set 
of physical assumptions concerning the origin, evolution, and equilibrium structure 
of halos in CDM universes. Observations of gravitational lenses have the potential 
to distinguish between the TIS profile and singular ones like the NFW profile, as 
several observable properties of gravitational lenses will be strongly affected by the 
presence, or absence of a central cusp in the density profile. One example of an 
important observable that can distinguish between various density profiles is the 
parity of the number of images. Lenses with nonsingular density profiles, such as 
the TIS, obey the odd number theorem. The number of images of a given source is 
always odd, unless the source is extended and saddles a caustic (see |124j . p. 172). 
Lenses with singular profiles, like the singular isothermal sphere, the NFW profile, 
or the Moore profile, need not obey this theorem, even for point sources. Most 
observed multiple-image gravitational lenses have either 2 or 4 images, and this 
may argue against profiles with a central core |121j . There are, however, other 
possible explanations for the absence of a third or fifth image. That image tends to 
be very close to the optical axis, and might be hidden behind the lens itself. Also, 
it is usually highly demagnified, and might be too faint to be seen. 

In this section, we derive all the lensing properties of the TIS. We also compare 
the TIS with three other density profiles: The Navarro- Frenk- White (NFW) den- 
sity profile, the Singular Isothermal Sphere (SIS), and the Schwarzschild Lens 6 . To 
compare the lensing properties of these various lens models, we focus on one par- 
ticular cosmological model, the currently favored COBE- normalized ACDM model 
with S7o = 0.3, Ao = 0.7, and Hq = 70kms~ 1 Mpc~ 1 (this model is also cluster- 
normalized). 



5.2 THE DENSITY PROFILES 
5.2.1 The Radial Density Profiles 

We will compute the lensing properties of halos with the TIS density profile, as well- 
approximated by equation l|61|l. and compare them with the properties derived for 
three comparison profiles. The first one is the Navarro, Frenk, and White (NFW) 
density profile, 

P(r) = 7- ^7 — y 2 , (79) 

V r / r NFvW v/ r NFW l ) 

where /?nfw * s a characteristic density and r NFW is a characteristic radius. We will 
also consider the Singular Isothermal Sphere (SIS) density profile, 

"M-a^a- (80> 

where oy is the velocity dispersion and G is the gravitational constant. This model 
might not represent actual halos very well, but it is a well-studied profile that has 

6 A point mass 
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important theoretical value. Finally, for completeness, we will also consider the 
Schwarzschild lens, 

p(r) = M Sch 5 3 (r) (81) 

where Mg c h is the lens mass and 5 3 is the three-dimensional delta function. 

In order to compare the predictions for halo lensing for different halo density 
profiles, we must relate the parameters which define one profile to those which define 
another. For this purpose, we shall use the same approach as |152j . The virial radius 
V2oo of a halo located at redshift z is defined, as usual, as being the radius inside 
which the mean density is equal to 200 times the critical density p c {z) = 3H 2 (z) /SttG 
at that redshift [where H(z) is the Hubble parameter]. The mass M200 inside that 
radius is given by 

Msoo = **>"frt« . (82 ) 

When comparing the lensing properties of different density profiles, we will consider 
halos that are located at the same redshift z = zl (the lens redshift) and have the 
same value of r2oo- By definition, these halos will also have the same value of M^oo- 
By stretching the terminology, we will refer to M200 as "the mass of the halo." 
This point needs to be discussed. For the Schwarzschild lens, M200 is indeed the 
mass of the halo. The SIS density profile drops as r~ 2 at large r, and the mass 
therefore diverges unless we introduce a cutoff. The halo mass will then be equal to 
M200 only if the cutoff is chosen to be r2oo- The NFW density profile drops as r~ 3 , 
hence the total mass diverges logarithmically, and this profile also needs a cutoff. 
The TIS density profile drops asymptotically as r~ 2 , but the TIS model includes a 
cutoff. This cutoff is located a radius r t ~ 1. 2^200 , and the mass inside the cutoff is 
Mi = 1.168M2oo- I n an Y case, a rigorous definition of the halo mass would require an 
unambiguous determination of the boundary between the halo and the background 
matter (such determination exists only for the TIS model), as well as dealing with 
the fact that the assumption of spherical symmetry that enters in these models most 
likely breaks down for real halos at large enough radii. Treating M200 as the actual 
mass of the halo is the best compromise. 

For a given halo mass M20O) redshift z, and cosmological background model, the 
density profiles for the SIS and Schwarzschild lens are fully determined. For the 
SIS, we integrate equation (|%Tj|) between r = and r = r2oo, and get 

M 20 o = ^|^. (83) 



Combining equations (|8*2*)) and (|83|). we get 



For the Schwarzschild lens, Mg c h is simply given by M2oo- 
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For the NFW profile, the concentration parameter c must be specified in addition 
to the parameters M200 an d z. The value of c is not completely independent of the 
other parameters since there is a statistical expectation that c is correlated with M200 
and z. However, for any individual halo, c is not known a priori. To determine c for 
a given halo, we will use the typical value expected from the statistical correlation of 
c with halo mass and the redshift of observation of the halo, z Q b s = zl,, according to 
the formalism of |44j . Once the value of c is known, we can compute the parameters 
of the profile. The characteristic radius of the NFW profile is given by 

r NFW = r 20o/c , (85) 

and the characteristic density /? NFW is given by 

_ 200cV(z) 
Pnfw ~ 3 [ln(l + C )-c/(l + c)] • (86) 

The TIS halo is uniquely specified by the central density po and core radius r$. 
These parameters are functions of the mass M200 and redshift z, but in this case, 
z is not the redshift zl where the halo is located, but rather the collapse redshift 
z co ii where the halo formed, which can be larger than z^. We will assume that a 
TIS halo which formed at some z co \\ > zl did not evolve between z co \\ and z^. For 
a given M200 and 2f co u, the parameters po and ro are determined as follows. The 
central density po depends only on z co ii, 

po = 1.8 x lOV(^coii) , (87) 

while M t = 1.167M 20 o = 772.6p ^, or 

r = 1.51 x l(T 3 (M 2 oo/po) 1/3 • 



Fig. 1341 shows a comparison of the various profiles, for halos with the same values 
of r2oo and M200 (the TIS curve is for z co \\ = z£). 

The equations describing the lensing properties of the Schwarzschild lens, the 
singular isothermal sphere, and the NFW profile can be found in the literature. For 
brevity, we will not repeat these equations here. We refer the reader to |124j for the 
Schwarzschild lens (§8.1.2) and the singular isothermal sphere (§8.1.4), and to j2H] 
and H52| for the NFW profile. 



5.2.2 The Interior Mass Profile 

The projected surface density is given by 



E(0 = / p(r)dz , (89) 
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Figure 34: Radial density profiles, for 4 different halos with the same values of 
r2oo and M2oo- Solid curve: TIS; dotted curves: NFW profiles with concentration 
parameters c = 2 and 10 (as labeled); dashed curve: SIS. 

where £ is the projected distance from the center of the halo, and z = (r 2 — ^ 2 ) 1 / 2 . 
For the TIS, we substitute equation (joTj) in equation and get 



J TIS 



(0 = TTPOfo 



A 



B 



(90) 



(a 2 r§ + e) 1/2 {b 2 r 2 + t 2 ) 1/2 

|68l I106j . For spherically symmetric lenses, one important quantity is the interior 
mass M(£) inside a cylinder of projected radius £ centered around the center of the 
lens. This quantity is given by 

M(0 = 2vr / 5 £(0^£'. (91) 
Jo 

we substitute equation (|90|) in equation (|91j) . and get 

M TIS (C) = 2vr 2 P0 r 3 [A(a 2 + e/r 2 ) 1/2 - B(b 2 + e/r 2 ) l/2 - Aa + Bb] . (92) 

Fig. 1351 shows a comparison of M(£) for various halos with the same values of r2oo 
and M 2 oo- 




Figure 35: Interior mass profiles, for 5 different halos with the same values of r2oo an d 
M2oo- Solid curve: TIS; dotted curves: NFW profiles with concentration parameters 
c = 2 and 10 (as labeled); short-dashed curve: SIS; long-dashed curve: Schwarzschild 
lens. 



5.3 THE LENS EQUATION 

Fig. illustrates the lensing geometry. The quantities rj and £ are the position of 
the source on the source plane and the image on the image plane, respectively, a 
is the deflection angle, and Dl, D$, and Dls are the angular diameter distances 
between observer and lens, observer and source, and lens and source, respectively. 
The lens equation is 

Dq 

V = j^C- D LS a (93) 

f |124j . eq. [2.15b]). Notice that since the lens is axially symmetric, we can write the 
quantities rj, £, and a as scalars instead of 2-component vectors. We nondimension- 
alize the positions and deflection angle, as follows: 



(94) 
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Figure 36: The lensing geometry: the dots indicate the location of the observer, 
lensing galaxy, and source. £ and r] are the positions of the image and the source, 
respectively, and a is the deflection angle. The angular diameter distances Dl, D^s, 
and Ds are also indicated. 



ro 
h 

D s r 



(95) 



a = DL ° LS& . (96) 



The lens equation reduces to 

y = x — a(x) . (97) 



"(*) = ■ 08) 



The deflection angle is given by 

M(r x 

f |124j . eq. [8.3]), where £ cr i t is the critical surface density, given by 

c 2 D s 



J crit 



4ttGD l D ls 



(99) 



where c is the speed of light. We substitute equation into equation (|9*H|) . and 

get 

a TIS (x) = ^-^x + ^ ~ ^ + ^ ~ M + Bb \ ' (1 ° 0) 

where k c is the central convergence of the TIS, defined by 

_£(£ = 0) wo M B 



^crit 

(see also |2B|). 



a 6 



(101) 
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5.4 CRITICAL CURVES AND CAUSTICS 

The lensing properties of a halo depend on its radial density profile p(r), the redshift 
zl where it is located, the source redshift zs, and the source position r/. In the 
remainder of this section, we will simplify the problem by fixing the redshifts zl 
and zs- We are primarily interested here in lensing by cluster-scale halos. Most of 
these halos are located in the redshift range 0.2 < Zl < 0.6, and none are located 
at redshift zl > 0.6. |151j In what follows, we will assume that halos are located at 
zl = 0.5. We will also assume that the source is located at redshift zs = 3. This 
is not a a very constraining assumption, because the lensing properties vary weakly 
with the source redshift for zs S> 1. 

For the Schwarzschild lens, SIS, and NFW profile, the density profile p(r) is fully 
determined by the values of M200 and Zl- For the TIS, we also need to specify the 
collapse redshift z co \\. We will consider three particular values, z co \\ = 0.5, 1.0, and 
1.5. 



5.4.1 Solutions 



The determination of the critical curves is quite trivial for axially symmetric lenses. 
Tangential and radial critical curves are defined respectively by 



m(x t ) 



a(x t ) 





A 


x t 


' d(m/x)~ 




da 


dx 




.dx _ 



1, (102) 
= 1. (103) 



(H2H eq. [8.3]), where m{x) = M(r x)/irr$ 

^crit is the dimensionless interior mass. 
For the Schwarzschild lens and the SIS, the solutions are xt = 1, and there are no 
real solutions for x r . For the NFW and TIS models, equations ()102|) and (|103|) have 
to be solved numerically for xt and x r . The solutions for the TIS are plotted in 
Fig. 03 as functions of k c . Also plotted is the radial caustic radius y r , obtained 
by substituting the value of x r into equation (|97|). (The value of y r we obtain is 
actually negative, because the source and image are on opposite sides of the lens. 
The actual radius of the caustic circle, then, is the absolute value of y r .) Both xt 
and \y r \ increase rapidly with k c , while the value of x r levels off. Fig. 1381 shows the 
angular radii of the tangential and radial critical circles, 9 r = ^ r /Di, and 9% = ^t/Di, 
in arc seconds as functions of the mass of the lens. For spherically symmetric lenses, 
multiple images (and thus critical circles) are possible only if the central convergence 
At(0) = £(0)/£ cr it exceeds unity ( |124| . p. 236, theorem [e]). For the Schwarzschild 
lens, SIS, and NFW profile, the central convergence diverges, hence these profiles can 
always produce multiple images. But for the TIS, the central convergence k(0) = k c 
is finite, and multiple images can be produced only if k c > 1. This explains the 
sharp low- mass cutoff seen in Fig. ISElfor the TIS (solid curves). 
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Figure 37: Radii of the radial critical circle, x r , tangential critical circle, xt, and 
radial caustic, y r , versus k c , for the TIS. 

The value of 6% for the Schwarzschild lens is called the Einstein radius Be- It is 
often used to estimate the characteristic scale of image features caused by strong 
lensing (e.g. ring radius, radial location of arc, image separations) and to estimate 
the size of the region within which the mass responsible for that strong lensing must 
be concentrated. Since lensing halos are not actually point masses, however, the 
angular radius B r { ng of the actual Einstein ring which results if the source is located 
along the line of sight through the lens center will usually differ from the Einstein 
radius Be, assuming that the lens mass distribution is actually capable of producing 
a ring. As we see in Fig. I38( 9e significantly exceeds Bt for all profiles considered 
(TIS, NFW profile, and SIS) for all masses considered. Hence, a mass estimate 
based on assuming that the scale of image features is of order 9e will underestimate 
the actual mass of the lens, unless the lens happens to be a Schwarzschild lens. 

A source located behind the lens will produce multiple images if y < y r . The 
angular cross section for multiple imaging is therefore 



In Fig. 1391 we plot the ratio of the cross sections for the NFW and TIS profiles. At 
low masses, M < 5 x 10 15 M o for z coU = 0.5, M < 1.2 x 1O 15 M for z co n = 1.0, 
M < 3 x 10 14 M Q for z co \\ = 1.5, the ratios are less than unity, indicating that a 
distribution of lenses described by the NFW profile would be more likely to produce 




(104) 
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Figure 38: Angular radius of the tangential critical curve (top) and radial critical 
curve (bottom), versus lens mass, for the TIS (solid curves), the NFW profile (dotted 
curves), the SIS (short dashed curves) and the Schwarzschild lens (long dashed 
curves). For the TIS, the values of z co \\ are indicated. 



cases with multiple images than if the same distribution is described by the TIS 
model. This trend is reversed at higher masses, and a TIS is more likely to produce 
multiple images than a NFW profile of the same mass. 
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Figure 39: Ratio of the cross sections for multiple imaging by the TIS and the NFW 
profile, versus lens mass. The values of z co \\ for the TIS are indicated. 

5.4.2 Illustrative Example 

Using a simple ray-tracing algorithm, we computed the image(s) of a circular source 
of diameter Ay = 1, created by a TIS with central convergence k c = 4.015. The 
results are shown in Fig. f° r 6 different locations of the source, ranging from 
y = 8.0 to y = 0.0. For each case, the left panel shows the source and the caustic 
circle (y r = 5.640) on the source plane, and the right panel shows the images(s), 
the radial critical circle (x r = 3.334), and the tangential critical circle (xt = 9.783) 
on the image plane. At y = 8.0, only one image appears. At y = 5.4, the source 
overlaps the caustic, and a second, radially-oriented image appears on the radial 
critical circle. At y = 4.8, the source is entirely inside the caustic, and the second 
image splits in two images, located on opposite sides of the radial critical circle, 
forming with the original image a system of 3 aligned images. As the source moves 
toward y = 0, the central image moves toward x = and becomes significantly 
fainter, while the other images move toward the tangential critical circle and become 
bright, elongated arcs. At y = 0, the two arcs have merged to form an Einstein ring 
located on top of the tangential critical circle, while the central image, very faint, is 
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Figure 40: Images of a circular source. Each pair of panels shows the source plane 
in the left panel, with the caustic, and the image plane in the right panel, with the 
radial (inner) and tangential (outer) critical circles. The position y of the source on 
the source plane is indicated. We used k c = 4.015, and a source of diameter Ay = 1. 



still visible in the center. 

5.5 IMAGE SEPARATION AND MAGNIFICATION 

The locations of the images are computed by solving the lens equation IJ97J) . For the 
TIS, this equation must be solved numerically. In Fig. 1411 we plot the separation 
between the two outer images as a function of the source location. The plot only 
extends to y/y r = 1, since larger values of y only produce one image. The solid and 
dotted curves show the separations for the TIS and NFW profile, respectively, with 
various values of k c and k s = /?NFW r NFw/^crit- The separation is fairly insensitive 
to the source location, and stays within ~ 15% of the Einstein ring diameter Ax = 
2xt for all values of k c and n s considered. The dashed line in Fig. ^2 shows the 
separation for the SIS, which is independent of y. For all profiles considered, the 
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Figure 41: Separation Ax between the two outer images, in units of the tangential 
critical radius x t , versus source location y in units of the caustic radius y r . The 
solid curves, from top to bottom, corresponds to TIS with k c = 10, 5, 2.5, and 1.2, 
respectively. The dotted curves, from top to bottom, corresponds to NFW profiles 
with k s = 1.0, 1.0, 0.5, and 0.2, respectively. The dashed line corresponds to the 
SIS. Results for the Schwarzschild lens are not plotted. 



image separation is always of order the Einstein ring diameter, independently of the 
source location. This is particularly convenient for theoretical studies, when the 
actual source location can be ignored |93j . 

The magnification of an image located at position x on the lens plane is given 

by 
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da 
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( |124j . eq. [8.17]). We computed the magnification of the images produced by halos 
of masses 10 13 , 10 14 , 3 x 10 14 , and 10 15 M Q . Fig. Ii2l shows the total magnification. 
The dotted curves show the results for a NFW profile. As y (or rj) decreases, the 
magnification slowly increases, until the source reaches the radial caustic y = y r . 
At that moment, a second image, with infinite magnification appears on the radial 
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critical curve (for clarity, we truncated those infinite "spikes" in Fig. 142(1. As y keeps 
decreasing, that second image splits into two images, and the total magnification 
becomes finite again, until the source reaches y = 0, and an Einstein ring with 
infinite magnification appears on the tangential critical curve. Of course, these 
infinite magnifications are not physical, since they can only occur for point sources. 
The total magnification is always larger than unity, and always larger when 3 images 
are present. 

The solid curves in Fig. 1421 shows the results for the TIS. At low masses, there is 
no radial caustic (See Fig. l38|) . and only one image appears. Because of the presence 
of a flat density core, the magnification is nearly constant if the path of the rays goes 
near the center of the core. For instance, for the case M200 = 1O 13 M0 (top left panel 
of Fig. 142(1 . tq > lOkpc, hence, over the range of r] being plotted, we are way inside 
the core. As the mass increases, the magnification increases, until a radial caustic 
forms. This happens at M 20 o = 1-11 x 1O 15 M for z co \\ = 0.5, M 20 o = 3.5 x 1O 14 M 
for z co n = 1.0, and M 20 o = 1.2 x 1O 14 M for z coU = 1.5 according to Fig. EHl At 
this point, the TIS has the ability to form three images, and the results become 
qualitatively similar to the ones for the NFW profile. 

The short-dashed and long-dashed curves in Fig. l42l shows the results for the SIS 
and the Schwarzschild lens, respectively. Because of the absence of radial caustic, the 
magnification always varies smoothly with source position, and the only divergence 
occurs at rj = 0, where an Einstein ring forms. 

For a given mass, the Schwarzschild lens always produces a stronger magnifica- 
tion than a SIS or NFW profile, unless the source is very close to the radial caustic 
of the NFW profile, where a "spike" of infinite magnification forms. The magnifica- 
tions produced by the SIS and NFW profile tend to be similar. If the NFW profile 
produces only one image (that is, we are on the right hand side of the dotted spike 
in Fig. I42|) . the SIS produces a larger magnification than the NFW profile. But if 
the NFW profile produces 3 images, then the total magnification exceeds the one 
produced by the SIS. As for the TIS, at low masses, where only one image forms, the 
magnification is much smaller than for the other profiles. But at large masses, where 
multiple images can form, the magnification becomes comparable to the one for the 
other profiles, and can even exceed the magnification produced by the Schwarzschild 
lens. 

5.6 WEAK LENSING 

Weak lensing usually refers to the magnification and distortion of the image of 
a background source by a foreground lens. Unlike strong lensing, weak lensing 
normally does not produce multiple images of single sources. The detection of 
coherent distortion patterns in the sky has been used to constrain the mass of 
clusters. The first detections were reported by |H3 El H32( I147j . followed by 
many others (see [131 and references therein). More recently, the distortion pattern 
produced by individual galaxies has also been detected [201 CUB 1121 EH1 EM [Ml H05j . 




Figure 42: Total magnification ^ t ot versus source location r\ for lenses of masses 
10 13 - 1O 15 M , for the TIS (solid curves), the NFW profile (dotted curves), the SIS 
(short dashed curves) and the Schwarzschild lens (long dashed curves). For the TIS, 
the curves, from top to bottom, correspond to z co n = 1.5, 1.0, and 0.5, respectively. 



Following the approach of |152j . we use as measure of the distortion produced 
by a lens the average shear 7 inside a distance £ = r2oo from the lens center. In 
practice, this quantity would be evaluated by averaging the shear of all images 
observed inside r2oo, after having eliminated foreground sources. We estimate this 
quantity by integrating the shear over the projected area of the cluster. The average 
shear inside radius x is given by 



2 f x 2 f x 
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where k(x') = S(x')/S cr it is the convergence. After some algebra, we get 
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Figure 43: Ratio of the average shear inside radius r2oo for the NFW profile and the 
TIS profile, versus mass of the halo. The values of z co \\ for the TIS are indicated. 
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We evaluate this expression at r = r2oo, or equivalently x = ^ooAo 
It reduces to 

7Tis( r 20o) = 408.67— , 



%00 



24.2. 



(108) 



-'crit 



where we used equation (|ST|) to eliminate po- 

In Fig. we plot the ratio 7nfw/7tis versus mass. The ratios are larger at 
lower collapse redshifts z co \\, and also decrease with increasing mass. This reflects 
the fact that as the mass increases, the concentration parameter c of the NFW 
profile decreases, while the ratio r2oo/ r o f° r the TIS remains fixed at 24.2. This 
figure can be qualitatively compared with the top left panel of Fig. 3 in |152j . We 
reach essentially the same conclusion as these authors, namely that using the av- 
erage shear to estimate the mass of lensing halos can lead to considerable errors if 
the wrong density profile is assumed. For z co \\ = 0.5, we find 7nfw/7tis > 1 a t 
small mass, and therefore the mass of a TIS would be underestimated if the lens 
is incorrectly assumed to follow a NFW profile. At high mass, the true mass of a 
TIS would be overestimated. For z m \\ = 1.0 and 1.5, 7nfw/7tis < 1 a t ah masses 
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considered, down to 3 x 10 Mq, and therefore the true mass of a TIS would always 
be overestimated. 

5.7 DISCUSSION 

In the previous two sections, we derived the effects of strong and weak lensing, 
respectively. Here, "strong lensing" refers to case with multiple images, arcs, or 
rings, while "weak lensing" refers to the magnification and shear of single images. 
We can divide the observed cases of strong lensing in two groups. The first group 
contains the "arc second" cases: multiple-image systems with image separations of 
order arc seconds, or rings with radii of that order [HI]- 7 In most cases, the lens 
is a single, massive galaxy, with possibly some additional contribution from the 
environment in which this galaxy is located ( |145j : see, however, |113| ^. A classic 
example is Q0957+561, the first gravitational lens to be discovered. The second 
group contains the "arc minute" cases, in which the lens is an entire cluster of 
galaxies. These lenses produce mostly giant arcs, with radii in the range 15" — 60" 
(see Table 1 of |151| ). The most famous case is the cluster CL 0024+1654, which 
produces multiple arcs. 

We showed that for all profiles considered, the image separation is weakly de- 
pendent on the source location (Fig. 141 j) . when multiple images actually form. If 
we neglect this dependence, the image separation AO i=a 20 1 can be read off the 
top panels of Fig. |3E1 We see immediately that galaxy-size objects cannot produce 
arc-second separations if they are described by the TIS or the NFW profile (under 
the assumption described in §2 that c for NFW halos is the typical value for halos 
at z bs = zl)- The TIS does not produce multiple images unless the mass is of order 
1O 14 M or above. The NFW profile can produce multiple images for any mass, but 
a separation AO > 1" requires a mass of order 1O 13 M0. Even the SIS needs a mass 
in excess of 10 12 M Q to produce arc-second separations. At the cluster scale, all pro- 
files are capable of producing arc-minute separations. In this limit, for a given mass, 
the separation is larger for the Schwarzschild lens and smaller for the NFW profile 
(because of the small concentration parameter at large mass). The separations are 
comparable for the SIS and TIS for z co \\ = zl, but for z co ii > zl the separation tends 
to be larger for the TIS than the SIS. As we indicated in the introduction, the TIS 
and NFW profile are applicable to dwarf galaxies and clusters of galaxies, but might 
not be applicable to ordinary galaxy-scale objects because baryonic processes are 
neglected. Hence, the inability of the TIS and NFW profiles to produce arc seconds 
separations with galaxy-size lenses is not a concern. 

Let us now focus on the TIS, and perform a simple calculation to estimate the 
probability that a TIS can produce multiple images. To compute k c , we substitute 
equations (|82[). (|87|). and (|5S)) in equation (jlOlj) . For the particular cosmological 
model (ACDM) and redshifts (zg = 3, zl = 0.5) we have considered in this paper, 

7 See |http://cfa- www.harvard.edu/castles] 
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where M15 = M200/IO M . Strong lensing requires re c > 1. Setting k c = 1, 
•Zcoii = ^Li and solving for M200, we get M200 = 1-102 x 1O 15 M . The requirement 
for strong lensing then becomes 



M15 > Mi 5 , : 



where 



M15, 



< 1.102. 



(110) 



(111) 



Hence, in the worst case scenario, a halo that formed at the latest possible collapse 
redshift z co \\ = zl, the required mass is about 1O 15 M . That value drops quite 
rapidly with increasing z co \\, down to about 1.2 x 1O 14 M for z co \\ = 1.5 (see Fig. OS)) . 
For the cosmological model we consider, 1-<t density fluctuations collapsing at red- 
shift Zcon = (0.5, 1.0, 1.5) have masses of about M 15 = (2 x 10~ 3 , 4 x 10" 4 , 9 x 10~ 5 ). 
Such "typical" objects will not be capable of producing multiple images of a source 
at redshift zs = 3, since the resulting values k c = (0.122,0.127,0.120) are smaller 
than unity. This simply indicates that multiple images are not produced by typical 
objects, which is certainly consistent with the fact that fewer than 30 arc- minutes, 
multiple-image systems have been observed. It is remarkable that the value of k c for 
typical objects is nearly independent of z co \\. It just happens the the dependences of 
Pc(z co ii) and M15 in equation ([109)1 nearly cancel out when M15 is set to the typical 
value at that redshift. 

Increasing k c above unity would require an object about 500 times more massive 
than a typical object at the same redshift. Objects of this mass are rare but do exist. 
We can make a simple estimate of how atypical such a massive object is. Over most 
of the mass range of cosmological interest (from small galaxies to clusters of galaxies) 
the CDM power spectrum can be roughly approximated by a power law P(k) oc k n , 
where k is the wavenumber and n ~ — 2. The rms density fluctuation 5 Tms is then 
given by <5 rms ~ k 3 ^ 2 P 1 ^ 2 (k) oc /c 1 / 2 . At a given redshift, different values of the 
wavenumber k correspond to different mass scales M according to M oc k~ 3 . The 
relation between rms density fluctuation and mass scale at fixed epoch is therefore 
approximated by 

5 rms oc M~ 1/6 . (112) 

Increasing the mass by a factor of 500 therefore reduces 6 vms by a factor of 500 1 / 6 ~ 
2.8. Because of the reduction in <5 rms , a 1-cr fluctuation (<5 — (J rms ) at this higher 
mass will no longer collapse by the same redshift (it will collapse later), but a 2.8-cr 
fluctuation (5 = 2.85 rnis ) will. Such fluctuations are rare, but not vanishingly rare. 
In Gaussian statistics, the probability that a randomly located point in space is 
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inside a 2.8-cr density fluctuation (i.e. 5 > 2.85 rms ) is about 1/200. Hence, one of 
every 200 halos would be capable of producing multiple images (of course, whether 
any halo actually produces multiple images depends on the location of the sources) . 

To give a specific example of a cluster-mass lens with z co \\ > zl, consider the 
case of the cluster CL 0025+1654. 126 showed that the mass-model derived by 
|146j to explain their lensing data for CL 0024+1654 at z = 0.39 is very well-fitted 
by a TIS halo with po « O.O64/i 2 M pc -3 and r$ f» 20/i _1 kpc. This central density 
implies that the halo collapse redshift is z co \\ sa 2.5 (i.e. z m \\ » zl). The halo mass 
in then M 20 o = 661.6p ro = 4.84 x 10 14 M o (for h = 0.7), and equation (fT09|) gives 
k c = 4.8 ^> 1. It is therefore not surprising that this cluster produces strong lensing. 

5.8 SUMMARY 

We have derived the lensing properties of cosmological halos described by the TIS 
model. The solutions depend on the background cosmological model through the 
critical surface density £ cr it, which is a function of the cosmological parameters and 
the source and lens redshifts, and the TIS parameters po and ro, which are functions 
of the mass and collapse redshift of the halo, and the cosmological parameters. 
By expressing the surface density of the halo in units of E cr it and the distances 
in units of ro, all explicit dependences on the cosmological model disappear, and 
the solutions are entirely expressible in terms of two dimensionless parameters, the 
central convergence k c and the scaled position y of the source. We have computed 
solutions for the critical curves and caustics, the image separations, and the total 
magnification. The ability of the TIS to produce strong lensing (multiple images 
and rings) depends entirely on k c . If k c < 1, only one image can form. If k c > 1, 
either one or three images can form, depending on whether the source is located 
outside or inside the radial caustic. When three images are produced, the central 
one is usually very faint, being highly demagnified. The angular separation between 
the two outermost images depends strongly on k c , but only weakly on the source 
location. 

For comparison, we derived (or extracted from the literature) the lensing prop- 
erties of three comparison models: the NFW profile, the singular isothermal sphere, 
and the Schwarzschild lens. Unlike the TIS, all of these profiles have a central 
singularity, which allows them to produce multiple images at any mass, provided 
that the source is sufficiently aligned with the lens. In practice, image separations 
large enough to be resolved can be achieved by galactic-mass objects only for the 
Schwarzschild lens, and by supergalactic-mass objects for all profiles. 

We applied the TIS model to the currently-favored ACDM universe, to calculate 
the central convergence k c expected for TIS halos of different masses and collapse 
epochs. We found that high-redshift sources (e.g. zs ~ 3) will be strongly lensed 
by TIS halos (i.e. k c > 1) only for cluster-mass halos, assuming that these halos 
formed at the redshift they are observed. As equation (|109|) shows, the halo mass 
required for strong lensing can be decreased by increasing the formation redshift 
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of the halo. However, as the formation redshift increases, the typical halo mass 
decreases, leading to a near-cancellation of the effect: the value of k c for a typical 
halo is about 0.12, independent of the collapse redshift. From this, we showed that 
a halo described by a TIS must form out of a ~ 2.8 — a fluctuation to be capable 
of producing strong lensing (k c > 1), no matter when the halo formed. Notice that 
weak lensing (magnification and shear of single images) does not require k c > 1, 
and can be caused by halos of any mass. 
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